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ABSTRACT 

* 

i 

! 

An  experimental  technique  is  utilized  in  which  a  liglTt-gas  » 

gun  is  used  to  launch  flat  impactor  plates.  to  high  velocities 
8  km/sec)  at  specimens  suspended  at  the  muzzle  of  the  gun. 
Impact-induced  shock  waves  at  pressures  to  ^  5  megabaro  are 
recorded  and  are  used  to  determine  the  shock  state  in  U  o 
specimen.  The  ability  to  launch  unshocked,  stress-fres  flat 
plates  over  a  wide  and  continuous  velocity  range,  coupled 
with  the  ability  to  launch  impactor  plates  of  the  same  material 
as  the  target,  results  in  hugoniot  measurements  of  relatively 
high  precision.  Measurements  were  made  on  Pansteel-77  (a  tung¬ 
sten  alloy),  aluminum  (2024-T4) ,  copper  (OFHC,  99.99%),  nickel 
(99.95%),  stainless  steel  (type  304),  titanium  (99.99%),  mag¬ 
nesium  (AZ31B) ,  beryllium  (S-200  and  1-400),  uranium  (depleted), 
Plexiglas,  and  quartz  phenolic.  The  results  are  compared  with 
those  of  other  researchers. 

Deviation  from  linear  shock  velocity  -  particle  velocity  was  1 

found  in  aluminum  beginning  at  v  1.0  megabars,  probably  attribu¬ 
table  to  melting  in  the  shock  front. 
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SECTION  I 

INTRODUCTION 


This  report  describes  the  experimental  procedures  and 
presents  the  results  of  the  work  performed  under  contract 
DA-18-C01-AMC-1126 (X) .  The  work  was  sponsored  by  the  U.  S. 

Army  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground, 
Aberdeen,  Maryland,  during  the  period  June  1966  to  June  1968. 
Very  high  pressure  measurements  were  made  to  determine  the 
hugoniot  equations  of  state  of  several  metals,  a  composite 
and  a  plastic.  Results  of  these  measurements  are  compared 
with  measurements  made  by  Al'tshuler,  McQueen,  Skidmore  and 
others.  This  work  extends  the  scope  of  a  General  Motors 
sponsored  research  project  reported  earlier ^  which  described 
the  use  of  a  light-gas  gun  to  obtain  pressures  substantially 
above  those  reported  by  other  researchers  in  this  country. 

The  materials  tested  are  copper  (OFHC,  99.99%),  nickel  (99.95%), 
titanium  (99.95%),  aluminum  (2024-T4) ,  stainless  steel  (type 
304) ,  magnesium  (AZ31B) ,  beryllium  (S-200) ,  beryllium  (1-400)  , 
Fansteel-77  (90%  W,  6%  Ni,  4%  Cu) ,  uranium  (depleted),  quartz 
phenolic,  and  Plexiglas  (Rohn  and  Haas  II  UVA) .  Table  1  lists 
the  measured  pressure  range  for  each  material. 

Section  II  of  this  report  presents  a  detailed  description  of 
the  experimental  techniques  employed  for  the  determination  of 
the  hugoniot  equations  of  state.  Data  analysis  is  briefly 
described  in  Section  III,  with  the  details  presented  in  Appen¬ 
dix  A.  Appendix  A  also  describes  and  lists  the  data  analysis 
computer  program  SHOVEL  used  to  reduce  the  experiment  data. 
Included  in  Appendix  B  is  a  summary  of  the  results  of  an  error 
analysis. 
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Experimental  results  and  comparison  with  other  experiments  are 
presented  in  Section  IV.  Section  V  contains  a  summary  of  the 
experimental  results  and  conclusions. 


TABLE  1 

SUMMARY  OF  MATERIALS  TESTED  AND 
PRESSURE  RANGES  EXAMINED 


PRESSURE  RANGE 

MATERIAL  _ (Mb) 


Fansteel-77  0.3  -  5.0 
Copper  (OFHC ,  99.99%)  1.0  -  4.5 
Aluminum  (2024-T4)  0.5  -  2.2 
Nickel  (99.95%)  0.8  -  4.7 
Stainless  Steel  (Type  304)  0.8  -  4.2 
Titanium  (99.95%)  0.4  -  2.7 
Beryllium  (1-400)  0.5 
Beryllium  (S-200)  0.3  -  1.6 
Quartz  Phenolic  1.3 
Magnesium  (AZ31B)  0.2  -  1.4 
Plexiglas  (Rohn  and  Haas,  II  UVA)  0.7  -  1.0 
Uranium  (depleted)  0.8  -  4.6 
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SECTION  II 

EQUATION  OF  STATE  MEASUREMENT 


In  the  past  two  decades,  high  explosives  have  been  used  to 

initiate  compressive  waves  with  amplitudes  from  tens  of  kilo- 

bars  to  several  megabars  in  many  materials.  High  explosive 

plane  wave  generators  placed  either  in  direct  contact  with 

the  material  or  in  contact  with  a  "standard"  material  upon 

which  the  specimens  were  placed,  were  used  to  measure  hugoniot 

states  to  approximately  600  kilobars.  Considerably  higher 

pressures  were  obtained  from  explosive  systems  in  which  the 

high  explosive  was  used  to  accelerate  a  thin  flier  plate 

across  a  gap  and  then  impact  the  specimen  surface.  In  this 

manner,  pressures  to  approximately  2  megabars  were  generated 

12) 

in  materials  of  high  density  by  McQueen  and  Marsh'  of 

(3) 

Los  Alamos  Scientific  Laboratory  (LASL) .  Hart  and  Skidmore 
increased  the  pressure  range  of  these  measurements  to  over 
5  megabars,  using  a  radially  converging  explosive  system  which 
accelerated  plates  to  high  velocities  at  some  expemse  in  pre¬ 
cision;  the  converging  shock  wave  system  adding  complexity  to 
.  (4) 

the  analysis,  Al* tshuler,  et  al,  extended  the  range  of 
measurements  to  above  10  megabars,  accelerating  his  flier  plates 
in  an  undisclosed  fashion.  Work  in  the  United  States  has  not 
progressed  above  the  2  megabars  reported  by  McQueen  until  this 
study,  which  provides  an  extension  of  Lower  pressure  data  to 
over  5  megabars. 

For  this  study  an  "accelerated  reservoir"  light-gas  gun  was 
used  to  accelerate  flier  plates  to  velocities  extending  above 
8  kilometers  per  second.  This  method  of  experimentation  offers 
significant  advantages  over  the  explosive  techniques  previously 
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used.  Unshocked,  stress-free  impactors  of  similar  or  dissimilar 
material  to  the  specimen  can  be  impacted  over  a  wide  and  con¬ 
tinuous  pressure  range.  Of  significance  is  the  simplicity  of 
the  calculation  of  the  shock  state  in  the  specimen,  using 
either  symmetric  impact  assumptions  or  the  measured  hugoniot 
of  the  gun- launched  impactor.  For  experiments  in  which  the 
shock  is  created  in  a  standard  by  either  direct  contact  with 
explosives  or  on  being  struck  by  an  explosively  accelerated 
plate,  the  hugoniot  point  is  less  readily  calculable.  In  this 
case,  it  is  necessary  to  assume  a  form  of  the  equation  of  state 
for  the  standard  in  order  to  get  the  hugoniot  point  for  the 
specimen. 


THEORETICAL  CONSIDERATIONS 


The  application  of  the  principles  of  conservation  of  mass, 
momentum,  and  energy  across  a  discontinuity  have  led  to  the 
well-known  Rankine-Hugoniot  equations.  The  equations  were 
derived  originally  for  fluids,  but  may  be  applied  to  solids 
when  the  pressure  P,  is  understood  to  represent  the  one-di¬ 
mensional  stress  normal  to  the  wave  front.  These  equations 
may  be  used  to  represent  the  discontinuous  change  of  pressure 
P,  density  p,  specific  volume  V,  and  internal  energy  E,  across 
a  shock  front  as  they  are  related  to  the  shock  wave  velocity 
Us,  and  the  particle  velocity  behind  the  shock  front  up, 
(Figure  1) . 


P 


p  U  u 
o  s  p 


(1) 


p  U  =  p(U  -u  ) 

OS  S  p 


(2) 


E-Eq  =  1/2  (P  +  Pq)  (VQ-V) 


(3) 
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Thus  a  measurement  of  the  shock  wave  velocity  and  the  particle 
velocity  associated  with  the  shock  wave  provides  sufficient  in 
formation  to  calculate  the  change  in  the  thermodynamic  state 
of  the  specimen  (assuming  that  steady  state  conditions  prevail 
behind  the  shock  front) . 


Material  At 
Shock  Conditions: 

Av  »*,,  o,,  e, 


Material  At 
Initial  Conditions: 


Shock 

Front 

Velocity  •  Us 

Figure  1  Schematic  of  Shock  Wave  Parameters 

Although  measurement  of  shock  velocity  is  relatively  straight¬ 
forward,  the  measurement  of  particle  velocity  is  more  diffi¬ 
cult  experimentally.  For  this  study,  two  techniques  were  used 
to  calculate  the  particle  velocity  associated  with  a  measured 
shock  velocity.  The  first  method  applies  by  symmetry.  For 
an  impact  of  a  specimen  launched  at  velocity  v  onto  a  specimen 
of  the  same  material  a  shock  wave  is  induced  with  particle 
velocity  equal  to  one-half  the  impact  velocity,  or 

Up  «  1/2  v  (4) 


5 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-68-13 

To  rigorously  apply  this  condition#  it  is  necessary  that  the 
impactor  and  specimen  be  in  the  same  thermodynamic  state# 
i.e.,  the  impactor  has  not  been  shock  heated  nor  subjected 
to  irreversible  changes  due  to  shock  loading  during  accelera- 
tion.  These  conditions  have  been  satisfied  for  this  study. 

Since  the  impact  velocity  is  measured  with  high  precision 
(customarily  ^  0.05%)#  the  particle  velocity  also  can  be 
calculated  with  similar  precision.  A  series  of  tests  are 
conducted  over  a  range  of  different  impact  velocities#  the 
highest  pressure  being  obtained  at  the  highest  impact  velo¬ 
city  (about  8  km/sec) .  Each  test  furnishes  a  point  on  the 
locus  of  final  compressed  states  known  as  the  hugoniot. 

Figure  2  shows,  in  the  pressure-particle  velocity  plane#  a 
graphical  description  of  the  conditions  of  symmetrical  impact. 


Symmetrical  Impedance  Match 
Schematic 
( Similar  Materials ) 


Figure  2  Schematic  of  Symmetrical  Impact  Analysis 
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To  obtain  pressures  higher  than  those  created  by  symmetrical 
impact  at  the  maximum  launch  velocity,  it  is  possible  to  im¬ 
pact  the  specimen  with  a  material  of  higher  shock  impedance 
(defined  as  the  product  of  the  initial  density  and  the  shock 
velocity)  and  to  calculate  the  resultant  particle  velocity  and 
pressure  by  an  adaptation  of  the  impedance  matching  technique 
developed  by  Walsh  et  al. ^  With  this  technique,  a  measure¬ 
ment  of  the  velocity  of  the  impactor  material,  the  hugoniot 
of  which  has  been  previously  measured,  is  sufficient,  when 
combined  with  a  measurement  of  the  shock  velocity  in  the  speci¬ 
men  material,  to  determine  a  point  on  the  hugoniot  of  the 
specimen.  Figure  3  shows  an  impactor  of  known  hugoniot  strik¬ 
ing  a  specimen  with  velocity  v.  A  single  shock  wave  of  velo¬ 
city  Us  is  induced  in  the  specimen.  The  intersection  of  the 
line  pQUg  and  the  hugoniot  of  the  impactor,  centered  at  velo¬ 
city  v,  determines  the  shock  pressure  and  the  particle  velocity 
in  the  specimen. 


Particle  Velocity 


Figure  3  Schematic  of  Dissimilar  Materials 
Impact  Analysis 
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EXPERIMENTAL  TECHNIQUES 

The  experimental  determination  of  hugoniot  equations  of  state 
using  the  impedance  match  technique  is  based  on  the  measure¬ 
ment  of  the  shock  velocity  in  the  specimen  and  of  the  velocity 
of  the  impactor.  With  the  techniques  described  below,  these 
two  measurements  can  be  made  with  precisions  of  approximately 
0.5  and  0.05  percent  respectively,  resulting  in  hugoniots  of 
good  accuracy  considering  the  limited  number  of  tests  conducted 
on  several  of  the  materials. 

The  light-gas  gun  range  and  basic  instrumentation  have  been 

(1  6  7  8) 

described  in  several  other  papers  '  '  '  '  '  and  are  again  in¬ 
cluded  here  for  completeness  of  this  report. 


Launching  Techniques 

The  gun  used  to  accelerate  the  impactor  is  an  accelerated- 
reservoir  light-gas  gun  with  a  launch  tube  bore  diameter  of 
either  29  mm  or  64  mm.  This  type  of  gun  maintains  a  reason¬ 
ably  constant  pressure  on  the  base  of  the  projectile  during 
the  launch,  allowing  a  relatively  gentle  acceleration  of  the 
impactor  materials. 

Figure  4  shows  the  layout  of  the  range.  The  gun  consists  of 
the  following  major  components: 

1.  Powder  chamber 

2.  Pump  tube,  89  mm  internal  diameter  by  12  m  long 

3.  Accelerated-reservoir  high-pressure  coupling 

4.  Launch  tube,  either  29  mm  internal  diameter  by 
8  m  long  or  64  mm  diameter  by  8  m  long 

5.  Instrumented  target  chamber  and  flight  range. 
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Figure  4  Layout  of  ARLG  Gun  Range 


When  the  gun  is  loaded  for  firing,  gunpowder  is  placed  in  the 
powder  chamber  and  the  pump  tube  is  filled  with  hydrogen.  The 
hydrogen  is  compressed  by  a  plastic  nosed  piston  which  has  been 
accelerated  by  the  burnt  gunpowder.  In  turn,  the  projectile  is 
accelerated  by  the  release  of  the  compressed  hydrogen  through 
a  high  pressure  burst  diaphragm. 


Prior  to  firing,  the  flight  range  and  instrument  chamber  are 

.  -2 

evacuated  and  then  flushed  with  helium  to  approximately  10 
Torr  to  eliminate  any  spurious  effects  due  to  gas  build-up  and 
ionization  between  projectile  and  target.  The  sealing  lips  on 
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the  rear  of  the  plastic  sabot  are  pressed  tightly  against  the 
sides  of  the  launch  tube  by  the  high  pressure  gas  and  effectively 
eliminate  blow-by  of  the  hydrogen  gas. 

Careful  attention  to  the  condition  of  the  launch  tube  is  nec¬ 
essary  for  successful  firing  in  the  high  velocity  ranges.  Bore 
linearity  of  better  than  0.2  mm  over  the  full  8  m  length  of  the 
launch  tube  is  maintained.  Internal  diameter  is  maintained  con¬ 
stant  to  within  0.01  mm.  Launch  tubes  are  cleaned  and  honed 
after  each  firing  and  are  removed  every  15  to  20  firings  for 
reconditioning . 

Figures  5  and  6  show  the  instrumentation  chamber  designed  for 
the  high  pressure  studies.  This  chamber  is  connected  to  the 
barrel  of  the  gun  through  an  O-ring  seal  to  allow  free  axial 
movement  of  the  launch  tube.  The  target  chamber  and  target  am 
shock-mounted  to  prevent  premature  motion  before  projectile  im¬ 
pact.  To  facilitate  this,  several  stages  of  mechanical  iso¬ 
lation  have  been  arranged  in  the  barrel,  I-beam  support  struc¬ 
ture  and  the  concrete  foundation. 

The  impact  chamber  is  a  steel  cylinder  of  61  cm  O.D.  and  1.5  m 
length.  Physical  access  and  instrument  ports  are  precisely 
machined  in  a  horizontal  plane  and  in  planes  45°  above  the  hori¬ 
zontal.  Two  stations  of  six  ports  each  are  accurately  spaced 
30.5  cm  apart. 

Operationally,  the  ports  are  closed  against  O-ring  vacuum  seals 
with  Plexiglas  or  magnesium  windows  for  optical  and  x-ray  access 
or  with  steel  cover  plates. 
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r.gujte  5  forget  Chamber  Set-up  with  Continous 
Wrining  Streak  Camera  in  Position 


STATION 

#1  STATION 


Figure  6 


Top  View  of  Target  Chamber 
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INSTRUMENTATION 

The  irapactor  velocity  measuring  system  consists  of  a  laser  trig- 

{ 8 ) 

gering  system  and  two  short  duration  ‘lash  x-rays.  With  this 
system,  impactor  velocities  are  measured  accurate  to  0.05%.  The 
triggering  system  consists  of  a  neon-helium  gas  laser  aimed  at 
a  photo-detector  across  the  impact  chamber  orthogonal  to  and 
intersecting  the  line  of  flight  of  the  projectile,  A  photo¬ 
multiplier  monitors  the  laser  light  output  through  a  set  of 
masks  and  a  narrow  band  optical  filter.  When  light  interruption 
occurs  due  to  projectile  passage,  a  sharp  change  of  voltage  level 
is  converted  into  a  signal  of  sufficient  amplitude  to  trigger  a 
Field  Emission  Corporation  30  nanosecond  dual  flash  x-ray  unit. 
The  x-ray  flash  exposes  a  Polaroid  film  plate  on  the  opposite 
side  of  the  chamber  by  means  of  a  fluorescing  intensifier  screen. 
The  trigger  and  x-ray  flash  system  is  tnen  duplicated  to  record 
the  passage  of  the  projectile  in  the  second  field  of  view  30.5  cm 
further  down  range. 


4 


The  spacing  between  the  two  x-ray  field  centerlines  is  indicated 
by  fiducial  wires  which  are  measured  by  an  optical  comparator  to 
within  0.2  iron.  Measurements  of  the  impactor  face  position  rela¬ 
tive  to  the  window  fiducials  allow  calculation  of  actual  pro¬ 
jectile  position  and  travel  over  the  time  interval  measured  be¬ 
tween  flash  exposures.  Figure  7  is  an  example  of  the  shadow¬ 
graphs  of  the  two  x-ray  stations  showing  the  projectile  in  free 
flight  before  impact. 


A  second  method  is  also  employed  to  measure  impactor  velocity. 

The  time  interval  between  the  first  x-ray  flash  and  the  impact 
of  the  projectile  on  the  target  is  recorded  electronically.  The 
impactor  and  target  positions  are  measured  from  the  x-ray  shadow¬ 
graphs  and  a  velocity  is  calculated.  Variations  in  measurements 
between  the  two  techniques  are  usually  less  than  0.05%. 
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Figure  7  X-ray  Shadowgraphs  of  Projectile 
Before  Impact 


The  target  is  located  approximately  60  cm  from  the  launch  tube 
muzzle  and  is  included  in  the  #2  x-ray  field  of  view.  Measure¬ 
ments  of  the  shock  wave  transit  time  in  the  target  are  made  us¬ 
ing  four  coaxial  self-shorting  pins  as  sensors .  The  shorting 
of  a  pin  results  in  a  sharply  rising  current  to  ground  which 
produces  a  signal  across  the  time-interval-meter  input  termina¬ 
tion  resistors.  The  circuit  is  so  designed  that  each  pin  signal 
can  be  seen  on  three  output  lines  and  is  free  of-  any  reflections 
or  ringing  for  several  hundred  nanoseconds.  The  individual  cir¬ 
cuits  are  "tuned"  by  the  use  of  trimmer  capacitors  so  that  the 
rise  time  of  each  signal  is  1.0  ±  0.1  nsec  to  12  volts.  Thus 
it  is  possible  for  the  combined  mechanical-electronic  signal 
system  to  make  use  of  the  ±  1/2  nsec  resolution  of  the  time  re¬ 
cording  instruments. 

The  shock  wave  transit  time-interval-meters  are  Eldorado  Model 
793  counters.  These  counters  have  a  specified  time  resolution 
of  ±  1/2  nsec  and  may  be  read  digitally  to  the  nearest  nano- 
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second.  They  require  an  input  signal  of  1  volt  with  a  rise 

time  approximately  1  nsec.  Although  instrument  stability  is 

4 

specified  to  be  one  part  in  10  for  long  term  and  five  parts 
in  10^  for  short  term,  in  actual  practice,  the  instruments  are 
calibrated  prior  to  each  shot  over  a  period  of  about  10  minutes 
The  shot  is  then  fired  within  five  minutes  of  completion  of 
the  calibration  procedure. 

The  planarity  of  the  shock  wave  induced  in  the  target  is  de¬ 
pendent  on  the  impactor  flatness  at:  impact.  The  impactor  sur¬ 
face  is  machine  lapped  and  then  hand  polished  flat  to 
0.5  x  10  3  mm.  Tests  performed  with  impactors  of  Fansteel-77 
and  OFHC  copper  indicate  the  surface  curvature  after  launch  to 
be  less  than  the  5  nanosecond  time  resolution  of  the  optical 
recording  system  at  a  launch  velocity  of  7  km/sec. 

The  impactor  tilt  relative  to  the  target  specimen  front  surface 
is  sensitive  to  launch  tube  linearity  and  sabot  alignment  as 
well  as  to  target  alignment.  The  capability  to  adjust  the 
target  position  and  perpendicularity  relative  to  the  launch 
tube  centerlii  a  by  an  optical  technique  brings  the  average  tilt 
at  impact  to  approximately  0.005  radians  (approximately  15  nano 
seconds  of  tilt  at  an  impact  velocity  of  7  km/sec) . 

Because  of  the  comparatively  gentle  acceleration  of  the  pro- 
jecti.e  to  its  terminal  velocity,  the  impactor  plate  is  not 
shock  heated.  In  addition,  free  flight  in  an  evacuated  range 
precludes  aerodynamic  heating.  This  accounts  for  the  flatness 
of  the  impactor  after  launch  and  significantly  reduces  the 
complexity  of  the  experiment.  The  estimated  temperatures  rise 
during  launch  of  the  order  of  1°C. 
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TARGET  DESIGN  AND  CONSTRUCTION 

The  29  mm  diameter  of  the  launch  tube  places  restrictions  on 
the  diameter  of  the  impactor  plate  and  on  the  size  of  the  speci¬ 
men  which  are  severe  enough  to  require  a  thorough  study  of  the 
optimization  of  target  dimensions.  In  general,  it  is  desirable 
to  operate  over  as  long  a  time  base  as  possible  for  transit  time 
measurements.  However,  the  launch  tube  diameter  controls  the 
allowable  specimen  thickness .  For  larger  specimens  rarefactions 
from  both  the  unconfined  edges  and  free  rear  surface  of  the  im¬ 
pactor  plate  could  overtake  the  head  of  the  shock  wave  before 
measurements  have  been  taken.  To  determire  the  maximum  speci¬ 
men  thickness  which  would  still  maintain  an  unrarefacted  area 
on  the  rear  surface  of  the  specimen  on  which  sensors  could  be 
placed  to  record  wave  arrival,  the  following  analysis  was  used. 

A  typical  estimate  of  the  angle  of  intrusion  of  plastic  rare¬ 
faction  waves  from  the  specimen  edges,  oi^,  is  to  assume  that 
tan  =  1  (see  Figure  8) ,  and  to  ignore  the  elastic  rarefac¬ 
tion  wave  system.  Although  for  many  materials  this  assumption 
is  justified,  for  some  materials  this  criterion  is  inadequate; 
in  particular,  materials  with  a  low  Poisson's  ratio  should  be 
calculated  more  carefully. 


The  elastic  wave  velocity,  Cg,  is  given  by 


C  =  C 
e  p 


^1=  C  K 

1+v  p 


where  v  is  Poisson's  ratio  and  is  the  plastic  wave  velocity, 
which  for  strong  shocks  is  a  function  of  u^,  the  particle  velo¬ 
city.  Existing  experimental  results  indicate  that  v  is  a  weak 
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Figure  8  Schematic  of  Target  and  Iirpactor 


function  of  the  shock  strength . To  a  good  approxima¬ 
tion,  the  elastic  rarefaction  angle,  a 2,  is 


tan  «2  =  K2tan2a^  +  (K2-l) 


/U  -u  \‘ 

M 


Figure  9  is  a  plot  of  a2  versus  Poisson's  ratio  for  the  metals 
listed  in  Table  2.  The  calculation  assumes  a  value  of  tan 


=  0.7,  taken  from  the  work  of  Al'tshuler 


,  who  notes  that 


for  compressions  greater  than  ^  1.3,  tan  becomes  essentially 
constant  at  0.70  ±  .03  and  Us  up  is  taken  as  unity  -  its  maxi¬ 


mum  vaJue.  Included  in  Figure  9  is  a  comparison  of  valves  of 
tan  a2  measured  in  this  laboratory  for  copper,  aluminum, 
titanium  and  beryllium.  As  these  measurements  fall  below  the 
calculated  line  of  the  minimum  allowable  design  angle,  it  is 
felt  that  equation  (6)  provides  a  reasonably  conservative  design 
critera.  In  practice,  the  design  angle  is  chosen  several  de¬ 
grees  larger  than  the  values  listed. 


16 


* w 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-68-13 


Figure  9  Intrusion  Angle  of  Elastic  Edge  Rarefaction 
vs  Poisson's  Ratio 

Impactor  thicknesses  were  chosen  to  avoid  rarefactions  originat¬ 
ing  at  the  impactor  free  rear  surface  from  overtaking  the  shock 
wave  until  shock  transit  time  measurements  were  complete.  A 
single  impactor  thickness  was  calculated  which  was  adequate  for 
all  materials  and  was  used  in  all  tests. 

The  target  specimen  was  a  machined  and  ground  disc  with  the  im¬ 
pact  and  rear  surfaces  machine  lapped  and  hand  polished  to  a 
surface  finish  of  1  microinch  rms  or  better.  The  lapping  pro- 

_3 

cedure  employed  produced  surface  flatness  within  ^  10  mm. 

-3 

Parallelism  was  maintained  to  'v  10  radians.  The  thicknesses 

-3 

of  all  specimens  were  measured  to  an  accuracy  of  ±  0.5  x  10  mm. 

The  target  specimen  thicknesses  at  the  pin  stations  were  mea¬ 
sured  with  a  Zeiss  light-section  microscope  employed  as  a  com- 
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para tor.  The  use  of  the  light-section  microscope  avoids  the 
problem  of  an  indicator  marring  the  specimen  surface  since  no 
physical  contact  is  made  with  the  surface.  Rather t  the  vertical 
position  of  a  thin  beam  of  light  projected  on  the  specimen  is 
compared  with  the  position  of  the  beam  projected  on  a  labora¬ 
tory  grade  gage  block  and  the  specimen  thickness  is  calculated. 

The  basic  features  of  the  target  design  employed  in  this  work 
are  illustrated  in  Figure  10.  Two  coaxial  shorting  pins  were 
passed  by  the  edge  of  the  specimen  disk  with  their  cap  faces 
exactly  in  the  plane  of  the  specimen  impact  surface.  These 
pins  were  used  to  initiate  the  timing  for  the  shock  wave  tran¬ 
sit  time  measurement  and  to  measure  impactor  tilt  in  terms 
of  the  time  interval  difference  between  their  respective  clo¬ 
sures  . 


Figure  10  Photograph  of  Target 
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TABLE  2 

CALCULATED  MINIMUM  DESIGN  ANGLES 
OR  MAXIMUM  EDGE  RAREFACTION  ANGLES 


CALCULATED  FROM  EQUATION  (6) 


Material 

Poisson's 

Ratio 

Minimum 
Design 
Angle,  a0 

Aluminum 

.332 

47° 

Beryllium 

.055 

o 

00 

IT) 

Copper 

.356 

45° 

Tantalum 

.342 

46° 

Tungsten 

.280 

50° 

Uranium 

.402 

43° 

Titanium 

.304 

O 

00 

Lead 

.430 

41° 

Magnesium 

.306 

49° 

Nickel 

.300 

00 

o 

Steel  (Mild) 

.290 

49° 

Plexiglas 

.327 

47° 

Two  rear  surface  pins  (or  one,  depending  upon  the  target  dia¬ 
meter)  were  mounted  in  line  with  the  tilt  pins  to  record  the 
shock  wave  arrival  at  the  rear  surface.  All  four  pins  were 
mounted  in  a  guide  fixture  which  assured  the  proper  geometrical 
spacing.  The  tilt  pins  were  fixed  in  position  in  a  dimension- 
ally  stable  epoxy,  while  the  rear  surface  pins  were  spring 
loaded  in  place  in  the  pin  guide.  The  pin  retainer  and  cable 
bracket  lent  rigidity  to  the  assembly  to  prevent  accidental 
damage  to  the  pin  shafts . 
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The  four-pin  targets,  in  conjunction  with  the  four  Eldorado 
one-nanosecond  time  interval  meters,  produced  four  values  of 
the  shock  wave  transit  time.  From  these  four  values  a  sin'-1's 
shock  wave  velocity  in  the  specimen  was  calculated  and,  by  a 
system  of  cross-checking,  an  indication  of  the  precision  of 
the  measurement  was  available. 

The  degree  of  non-planarity  of  impact  between  target  and  impactor 
is  calculated  by  comparing  shock  transit  times  recorded  by  the 
counters  started  by  each  of  the  two  front  surface  pins.  The 
time  difference  between  the  shorting  of  the  front  surface  pins, 
Atfc  (which,  when  combined  with  the  impact  velocity,  yields  the 
tilt  angle,  0)  is  calculated  from  At,  =  t.  -t„  _  =  t.  -,-t,. 
where  t^_c  is  the  time  recorded  on  the  counter  started  by  pin  #1 
and  stopped  by  pin  C  (Figure  10a) . 

For  the  highest  velocity  tests  (7-8  km/sec) ,  it  was  necessary 
to  lighten  the  projectile  by  reducing  the  diameter  and  the  thick¬ 
ness  of  the  impactor  plate.  The  target  designed  for  these  high¬ 
est  pressure  shots  had  a  slightly  smaller  diameter  and  thickness 
and  was  provided  with  only  one  coaxial  shorting  pin  on  the  rear 
surface . 

A 

The  coaxial  self-shorting  pins  employed  in  this  work  as  sen¬ 
sors  consisted  of  a  one  millimeter  diameter  tube  of  brass  sur¬ 
rounding  a  teflon  sleeve  and  a  copper  inner  conductor.  The 
pins  were  connected  to  RG174  50ft  cable  by  soldered  joints  and 
were  made  self-shorting  by  the  placement  of  a  brass  cap  over 
the  sensing  end,  which  left  a  small  gap  (on  the  order  of  0.050 
±  .002  mm)  between  its  inner  face  and  the  flat  end  of  the  inner 
conductor.  When  a  large  amplitude  stress  wave  reaches  the  cap 
face,  the  cap  is  set  into  motion  and  the  gap  is  closed  at  the 

A 

Model  CA-1039,  Edgerton,  Germeshausen  and  Grier,  Santa 

Barbara,  California. 
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free  surface  velocity  of  the  cap  material.  The  pin  gaps  were 
measured  by  x-ray  shadowgraphy ,  of  which  Figure  11  is  an  example 
and  the  measurements  were  employed  in  the  shock  velocity  calcu¬ 
lations  for  corrections  to  pin  closure  times. 


Figure  11 


X-ray  of  Coaxial  Shorting  Pin 
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SECTION  III 

DATA  ANALYSIS 


The  analysis  of  the  experimental  data  obtain  .  \n  the  research 
program  is  base  '  upon  the  impedance  match  solution  for  the 
determination  of  particle  velocity,  pressure,  energy  and  volume 
of  the  shock  state  in  the  specimen.  The  requisite  information 
in  the  analysis  is  the  shock  wave  velocity,  the  impact  velocity, 
and  the  hugoniot  of  the  impactor.  A  general  description  of  the 
method  of  data  analysis  is  given  here,  with  the  detailed  equa¬ 
tions  presented  in  Appendix  A. 


SHOCK  WAVE  VELOCITY 

The  measurements  relevant  to  shock  wave  velocity  are  the  target 
thickness  and  the  shock  wave  transit  time  interval.  In  order  to 
calculate  the  shock  wave  transit  time,  it  is  necessary  to  make 
refinements  upon  the  recorded  time  interval. 

Two  sources  of  refinement  to  the  measurement  are: 

(i)  Inclusion  of  the  effects  of  shock  wave  tilt  resulting 
from  rrn-planar  projectile  impace  on  the  target. 

(ii)  Correction  for  the  differences  of  closing  times  of 
coaxial  pins  with  different  gaps,  which  involves: 

(a)  Calculation  of  the  interaction  of  the 
impactor  with  the  two  front  surface  pins 

(b)  Calculation  of  the  interaction  of  the 
specimen  material  with  the  two  rear  sur¬ 
face  pins. 
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To  calculate  the  closing  velocity  of  the  cap  for  the  direct 
impact  of  the  projectile  material  on  the  front  surface  pins, 
an  impedance  match  solution  is  applied,  using  the  impact  velo¬ 
city,  the  hugoniot  of  the  impactor  and  the  hugoniou  of  the  cap 
material  (brass) . 

The  pin  gap  correction  for  the  interaction  of  the  rear  surface 
pins  and  the  specimen  is  based  on  the  impedance  match  solution 
of  the  shock  wave  in  the  target  being  transmitted  into  the  pin 
material.  The  hugoniot  of  the  specimen  must  first  be  estimated 
to  provide  the  necessary  constants.  The  shock  transit  time  is 
first  calculated  with  no  pin  gap  corrections  and  a  preliminary 
hugoniot  point  is  determined  for  the  specimen.  This  hugoniot 
point  is  then  used  to  calculate  pin  interactions  and,  through 
a  series  of  iterations  (usually  two)  the  preliminary  hugoniot 
point  is  modified  until  satisfactory  convergence  is  reached 
(differences  <  0.01%). 


IMPACT  VELOCITY 

Measurement  of  impact  velocity  has  been  discussed  earlier  under 
"Instrumentation".  The  flash  x-ray  system-used  furnishes  high 
precision  velocity  determination  providing  the  projectile  main¬ 
tains  a  constant  velocity  during  the  time  of  measurement.  A 
check  on  this  premise  is  provided  by  the  seconding  system  which 
measures  velocity  over  a  longer  baseline.  No  evidence  of  pro¬ 
jectile  acceleration  or  deceleration  during  its  free  flight 
has  been  observed. 

HUGONIOT  OF  THE  IMPACTOR 

Measurement  of  impactor  hugoniots  for  tests  in  which  specimens 
are  impacted  with  dissimiliar  materials  are  discussed  in  the 
following  section  "Experimental  Results". 
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SECTION  IV 
EXPERIMENTAL  RESULTS 


The  experimental  work  effort  was  first  concentrated  on  the 
measurement  of  the  hugoniots  of  the  three  materials  to  be  used 
as  impactor  plates  for  dissimilar  material  impacts.  The  ma¬ 
terials  we re  chosen  to  (i)  cover  a  range  of  shock  impedances, 
(ii)  be  easily  obtainable  and  consistent  in  their  material  pro¬ 
perties  from  batch  to  batch,  and  (iii)  to  coincide  with  stand¬ 
ards  chosen  by  other  laboratories .  The  materials  investigated 
were  Fansteel-77  (a  tungsten  alloy) ,  OFHC  copper,  and  2024-T4 
aluminum.  The  ratio  of  shock  impedances  of  the  copper  and 
Fansteel-77  with  respect  to  the  aluminum  is  approximately  1:2:4. 
These  three  materials  were  more  thoroughly  investigated  than 
the  remaining  materials  sc  that  errors  in  the  impedance  match 
solution  for  materials  impacted  by  these  standards  would  be 
minimized. 

A  typical  test  series  for  the  determination  of  the  hugoniot  of 
a  specimen  material,  for  instance  nickel,  began  with  a  series 
of  shots  using  nickel  for  both  impactor  and  target  over  the 
full  velocity  range  of  the  gun.  To  obtain  pressures  higher 
than  those  created  by  like-like  impact  at  the  highest  velocity, 
the  series  continued  with  impacts  using  a  material  of  higher 
impedance  than  the  specimen;  in  the  case  of  nickel,  Fansteel-77 
was  used.  Impact  velocities  were  adjusted  to  space  the  shots 
over  the  pressure  range  to  be  investigated. 

In  the  following  section,  the  experimental  data  are  presented 
in  tabulated  and  graphical  form.  Fits  to  the  shock  velocity 
vs  particle  velocity  data  have  been  made  by  the  method  of  least 
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squares  and  are  listed  for  each  material.  The  tabulations  in¬ 
clude  measured  and  calculated  parameters  and  an  indication  of 
the  "weighting  factor"  used  in  the  least  square  fits.  In 
general,  test  results  were  weighted  according  to  whether  the 
target  had  single  or  double  pins,  the  double  pin  targets  gen¬ 
erally  having  a  higher  weighing  factor  due  to  the  redundant 
measurements  of  shock  wave  velocity. 

The  tabulated  data  also  include  the  impactor  material  and  the 
impact  velocity.  Additional  figures  are  included  to  illustrate 
comparison  with  other  researchers . 

FANSTEEL-77 

Fansteel-77,  a  tungsten  alloy  composed  nominally  of  90%  tungs¬ 
ten,  6%  nickel,  and  4%  copper  was  employed  as  the  standard  for 
the  highest  pressure  tests.  Fansteel-77  was  chosen  over  pure 
tungsten  because  the  metallurgical  structure  reduces  the  brittle¬ 
ness  of  the  material  (which  can  cause  plate  fracture  during 
launch)  while  maintaining  high  strength  and  density. 

In  the  initial  work  the  quality  control  of  the  Fansteel-77 
stock  material  presented  several  problems.  The  material  is 
produced  by  powder  metallurgical  techniques  which  include 
pressing  and  sintering  of  a  billet  of  the  material.  The 
porosity  of  the  surface  of  the  billet  was  found  to  be  somewhat 
dependent  upon  the  compacting  pressure  prior  to  sintering. 

The  core  of  the  billet,  however,  was  found  upon  metallographic 
examination  to  exhibit  essentially  no  porosity,  and  sample  to 
sample  variations  in  density  were  less  than  0.5%  when  the  out¬ 
side  1  mm  was  removed  from  a  50  mm  diameter  bar.  The  chemical 
and  physical  properties  of  Fansteel-77  are  presented  in  Table  3. 
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TABLE  3 

CHEMICAL  AND  PHYSICAL  PROPERTIES 
OF  FANSTEEL-77 


Chemical  Properties 

Element 

W 

Cu 

Ni 


Wt  % 
90  ±  .1 
3.8  ±.6 
6.1  ±  .2 


Physical  Properties 


Yield  Strength  (0.2%  elong. ) 

Ultimate  Tensile  Strength 
Density 

Poisson's  Ratio: 

* 

Acoustic  Velocities  :  Longitudinal, 

Shear  , 

Bulk  , 


85,000  psi  (min) 
98,000  psi  (min) 

3 

17.01  ±  .01  gm/cm 
0.286 

C^  =  5.049  km/sec 

C  =  2.765  km/sec 

s 

CQ  =  3.912  km/sec 


Ultrasonic  tests  were  performed  by  J.  Havens  and 
R.  Lingle  of  this  laboratory. 


v 
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The  test  series  on  Fansteel-77  resulted  in  a  hugoniot  over  the 

pressure  range  .35  Mb  to  5.0  Mb.  Symmetric  impacts  were  used 

on  all  tests  so  that  particle  velocity  may  be  assumed  to  be 

one-half  the  impact  velocity.  Both  second  and  third  order  fits 

were  made  to  the  data;  however,  the  linear  relationship 

U  =  C  +  Su  describes  the  data  with  the  root  mean  square  (RMS) 
s  p 

deviation  not  significantly  larger  than  the  higher  order  fits. 
The  relationship  is  given  by: 

U  =  4.008  +  1.262  u  km/sec 
s  p 

RMS  deviation  of  U  was  ±  0.021  km  per  second  for  the  seventeen 

s 

data  points.  The  results  are  presented  in  Table  4. 

The  data  points  taken  are  displayed  in  the  Ug  vs  u^  plane 
in  Figure  12.  In  one  shot  the  velocity  was  not  measured 
directly  with  the  x-ray  system  and  instead  was  calculated 
from  the  gun  firing  parameters.  Estimates  made  in  this  man¬ 
ner  are  quoted  with  standard  errors  in  velocity  of  ±  2%. 

In  Figures  13,  14  and  15  are  displayed  the  data  for  Fan¬ 
steel-77  compared  with  the  data  obtained  by  researchers  of 

the  Ballistics  Research  Laboratory  of  the  United  States  Army, 

*  .13) 

Aberdeen  Proving  Grounds  and  the  data  of  Hart  and  Skidmore 

for  a  tungsten  alloy  similar  to  Fansteel-77. 

The  Ballistics  Research  Laboratories  data  are  in  excellent 
agreement  with  the  results  obtained  here.  The  hugoniot  mea¬ 
sured  by  Hart  ar 1  Skidmore  has  the  quadratic  form: 

U  =  2.95  +  2.47  u  -  0.342  u  2  (km/sec) 
s  P  P 

* 

Mr.  George  Hauver,  private  communication. 
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Shock  Velocity  vs  Particle  Velocity  for  Fansteel-77 
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Figure  13 


PRESSURE 
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MATERIAL  -  FANSTEEL 


SPECIFIC  VOLUME  -  CC/GM 


0-09 


Figure  15 


Pressure  vs  Specific  Volume  for  Fansteel-77 


MANUFACTURING  DEVELOPMENT  e  GENERAL  motors  corporation 


MSL-68-13 

and  falls  below  the  present  data  at  lower  values  of  u^.  This 
behavior  is  not  readily  explainable  since  the  intercept  of 
the  linear  fit  to  the  present  work,  CQ  =  4.008  km/seo.,  is  in 
agreement  with  the  bulk  sound  velocity  measured  ultrasonically , 
CQ  =  3.912  km/sec. 

It  is  interesting  to  note  that  the  hugoniots  for  pure  tungsten 
(as  measured  by  LASL)  and  for  Pansteel-77  are  quite  similar, 
the  LASL  linear  fit  for  tungsten  being  given  by; 

U  =  4.029  +  1.237  u  km/sec 
s  p 

although  the  density  of  the  tungsten  is  13%  higher  than  that  of 
Fans tee 1-7 7. 


OFHC  COPPER 

OFHC  copper  of  99.99%  purity  is  employed  as  a  standard  for 
the  intermediate  and  high  pressure  ranges  in  the  hugoniot 
experiments  of  a  number  of  laboratories.  Physical  character¬ 
istics  of  the  copper  are  described  in  Table  5  The  results  of 
twelve  tests  are  presented  here  in  Table  6.  The  linear  fit  to 
the  data  in  the  U  -  u  plane  is  given  by; 

®  hr 

U  =  3.964  +  1.463  u  km/sec 
s  p 

Root  Mean  Square  Deviation  of  Ug  =  i  0.009  km/sec  for  12  data 
points.  The  hugoniot  equation  recently  published  by  LASL^12^ 
is  given  by; 

U  =  3.940  +  1.489  u  km/sec 
s  p 

in  close  agreement  with  the  present  results. 
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TABLE  5 


CHEMICAL  &  PHYSICAL  PROPERTIES 
OF  OFHC  COPPER* 


Chemical  Composition 


Element 

Wt  % 

Iron 

0.0005 

Sulphur 

0.0025 

Silver 

0.0010 

Nickel 

0.0006 

Antimony 

0.0005 

Lead 

0.0006 

Copper 

Remainder 

Properties 


Density 

Poisson's  Ratio: 


8.930  gm/cm~ 
0.332 


Acoustic  Velocities:  Longitudinal,  C.^  =  4.757  km/sec 

Shear  ,  Cg  =  2.247  km/sec 

Bulk  ,  C  =  3.99  km/sec 


,  CQ  =  3.99  km/sec 


*  Based  on  manufacturers  specifications. 
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The  data  are  presented  in  Figure  16/  which  now  includes  a  legend 
to  indicate  the  impactor  material  used.  The  use  of  different 
impactors  provides  a  good  method  for  cross-checking  the  hugoniots 

(ID 

of  the  standards  as  suggested  by  McQueen  .  The  hugoniot  es¬ 
tablished  by  Fansteel-77  impacting  copper  is  indistinguishable 
from  the  copper  on  copper  tests.  This  both  checks  the  accuracy 
with  which  the  Fansteel-77  hugoniot  was  determined  and  provides 
additional  confidence  in  the  use  of  impedance  matching  tech¬ 
niques  as  used  in  this  experimental  system.  A  single  test  was 
conducted  using  2024  aluminum  as  an  impactor.  All  tests  were 
used  in  the  calculation  of  the  fit  given  above.  The  resulting 
hugoniot  displays  a  very  small  RMS  deviation/  0.1%  in  shock 
velocity) . 

-  .  \ 

Figures  17,  18  and  19  compare  the  present  data  with  the  data 
of  A1 '  tshuler  ^ ,  Walsh  ^  and  McQueen^.  Agreement  within 
'v  2%  in  the  linear  fit  was  obtained  (i.e.,  the  reported  shock 
velocities  are  within  'v  2%  of  the  values  predicted  by  the  pre¬ 
sent  linear  fit) . 

A  close  examination  of  the  data  indicates  a  small  (<  1%)  devi¬ 
ation  from  a  linear  fit  beginning  at  u  =  2.4  km/sec.  It  is 

A 

likely  that  this  deviation  is  associated  with  melting  in  the 
s  sk" 'front  (see  the  following  discussion  or.  2024-T4  aluminum). 
This  phenomenon  will  be  discussed  more  fully  in  a  forthcoming 
report. 

2024-T4  ALUMINUM 

The  hugoniot  e cperiments  performed  on  2024*^4  aluminum  extend 
over  a  pressure  range  of  0.45  to  2.2  Mb.  Fourteen  tests  were 
conducted,  employing  the  following  impactors:  Fansteel-77  - 
5  tests;  OFHC  copper  -  5  tests;  and  2024-T4  aluminum  ■  4  tests. 

In  addition,  Shot  No.  98  from  the  series  on  copper  is  included 
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Shock  Velocity  vs  Particle  Velocity  for  Copper 
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Figure  18  Pressure  vs  Particle  Velocity  for  Copper 
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Figure  19  Pressure  vs  Specific  Volume  for  Copper 
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as  a  cross-check  by  reversing  the  impactor  and  specimen  materials 
and  impacting  an  aluminum  plate  into  a  copper  target.  If  the 
hugoniot  of  the  copper  is  assumed  to  be  known  the  state  in  the 
aluminum  may  be  calculated. 

The  chemical  and  physical  properties  of  2024-T4  aluminum  are 
tabulated  in  Table  7  and  the  data  are  presented  in  tabular  form 
in  Table  8^  The  measured  initial  density  of  the  aluminum  was 
2.783  ±  .001  gm/cm3* 

Figures  20  and  21  are  a  plot  in  the  Ug  -  up  plane  of  th* ?  ex¬ 
perimental  results.  The  linear  fit  is  from  LASL,  Reference  22. 

A  departure  from  linear  behavior  is  seen,  beginning  at 
U  —  3.5  km/sec  (P  *  1  megabar),  where  shock  velocity  falls 
below  the  line  representing  a  linear  fit  to  the  data.  Although 
the  data  shows  scatter,  it  is  felt  that  the  trend  of  the  data 
in  this  region  is  beyond  experimental  error.  The  linear  fit  of 
LASL,  U  =  5.328  +  1.338  u  ,  to  data  below  this  pressure  range, 
was  compared  to  the  present  data  by  calculating  deviations  of 
the;  data  from  this  fit.  The  deviations  are  plotted  in  Figure  22 
along  with  other  high  pressure  data  from  Russian  researchers. 
Only  the  ten  tests  showing  least  internal  scatte .  and  tilt  are 
plotted. 


Urlin^15^  has  proposed  a  model  for  melting  in  the  front  of  a 
shock  wave  and  predicts  an  observable  effect  on  the  linear 
Us  -  Up  relation.  For  aluminum  the  melting  is  calculated  to 
begin  at  approximately  1  megabar.  The  present  data  follows 
the  tregd  predicted  by  Urlin,  although  whether  melting,  ex¬ 
perimental  inaccuracy,  or  other  phenomena  is  the  explanation 
for  the  large  deviations  found  between  =  3.5  to  4.5  km/sec 
remains  to  be  verified. 

The  comparison  of  hugoniot  data  from  other  workers (9,13,14,16) 
is  shown  in  Figures  21,  2.3  and  24,  and  displays  the  area  of 
divergence . 
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If  all  the  present  data  points  are  combined  the  data  may 
be  represented  by  the  equation: 


U  =  5.471  +  1.310  u 
s  p 


Root  Mean  Square  Deviation  =  ±  .022  km/sec  for  11  data  Faints.. 


TABLE' 7 


CHEMICAL  6  PHYSICAL  PROPERTIES 
OF  2024-T4  ALUMINUM 


Chemical  Composition 


Element 


Wt  % 


Silicon 

Iron 

Copper 

Manganese 

Magnesium 

Chromium 

Zinc 

Aluminum 


0.5 

0.5 

3.8  -  4.9 
0.3  -  0.9 
1.2  -  1.8 
0*10 
0.25 

Remainder 


Physical  Properties 


Yield  Strength  47,000  psi 

Ultimate  Tensile  Strength  68,000  psi 

Hardness  (Brinell  No.)  120 

3 

Density  (measured)  2.783  gm/cm 

Poisson's  Ratio:  0.332 

Acoustic  Velocities:  Longitudinal,  C^  =  6.38  km/sec 

Shear  ,  C_  =  3.20  km/sec 


2.783  gm/cm* 
0.332 


,  C  =  3.20  km/sec 
s 


Bulk 


,  CQ  =  5.20  km/sec 


Based  on  Alcoa  Aluminum  Handbook  and  specimen  certification. 


II 


> 


asm 


Shot  Impactor  Impact,  ,  Shock  Particle  Pressure  Relative  Volume  Density  Weighting 
Material  Velocity'  Velocity  Velocity  Volume  Factor 

(km/sec)  (kt^sec)  .  (kn^sec)  (Mb)  (cia^/gm)  (gin/cm-1) 


C-1238  FS  2 . 583  8.099  2.017  0.455  .?510  .2698  3.706 
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PARTICLE  VELOCITY  -  KM/SEC 


RATION 


Figure  20 


Shock  Velocity  vs  Particle  Velocity  for 
2024-T4  Aluminum 
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Figure  21  Shock  Velocity  vs  Particle  Velocity  for 
Aluminum 
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MATERIAL  -  ALUMIMW 
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Figure  23  Pressure  vs  Particle  Velocity  for  Aluminum 
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Figure  24  Pressure  vs  Specific  Volume  for  Aluminum 
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DEPLETED  URANIUM 

Five  test's  were  performed  for  the  determination  of  the  hugo- 
rtiot  of  depleteduraniumover  the  pressure  range  0 .9  to  4.6 
megatbarsi  The  impactors  were  either  OFHC  copper  or  Fan- 
steeitM.  No-  similiar  material  impact  tests  were  performed.  . 
^The  cfa^f  is’  {shipwh  ih.  Figure s  25 ,  26  and  27 .  Excellent  agree- 
iteht^wis  obtained;  wi  th  the  high  pr  as  sure  (above  2  megabars) 
data  oA  ^Skxtimore  ;and rMori-is . 

~  «  •  —  "*5  t ~  «. 

The  density  ofthe  sampleswas  18 .951  gm/cm2 .  The  chemical 
ahd  phyi3ical:  prppbrties  are  summarized  in  Table  9,  and  the  hug- 
bnidt  data  are  presented  in  Table  10.  It  was  noted  that  the 

lapped  eurf aces  pXidi zed  rapidly ,  turning  from  a  light 
tiro  -€p  a  da^ybluetbfown  color.  That  this  oxide  layer  is  very 
thin  may  be  demonstrated  by  the  fact  that  several  swipes  of 
the  surface  on  awetiapp  removes  the  darkened  layer. 

The  1  in'e’ar  fit  to  the  Resent  dita  is  given  by  the  equation : 

U  =  2;  443  +  1.582  u-  km^sec 

s  ------  .  — -  -:p- -  - 

-Root  Mean  Square  Deviation-  ±  ;029  km/sec  for  5  data  points. 

Skidmore's  data  for  5  data  points  is  included  in  the  figures 
and  is  fit  by  the  .quadratic  equation : 

U  =  2.55  +  1.504  U  +  .0901  (u  -  7.50) 2  km/sec 
s  P  P 

for  u_  >  2.5  km/se 3.  ^  ~ 

tor  ~  ~~ 

The  linear  fit  from  LASL^2^  is  D-  =  2.487  +  1.539  u  . 

.  -  .  .  s  p 
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Figure  25 


Shock  Velocity  vs  Particle  Velocity  for 
Ura.dum 
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TABLE  9 

CHEMICAL  AND  PHYSICAL  PROPERTIES  OF  DEPLETED  URANIUM 


Chemical  Properties 


Element  Wt  (ppm) 

Aluminum  9 

Boron  0 . 3 

Cadmium  0 . 3 

Chromium  3 

Copper  8 

Iron  10 

Magnesium  2 

Manganese  15 

Molybdenum  2 

Nickel  20 

Lead  1 

Silicon  40 

Samarium  1 

Uranium  wt  %  99.8% 


Physical  Properties: 

Yield  (0.1%  elongation) 

Tensile  Strength 

Density  (measured 

Poisson's  Ratio: 
Acoustic  Velocities: 


47,250  psi 
124,000  psi 
18.951  gm/cm3 
0.402 

Longitudinal,  C^  =  2.97  km/s-sc 
Shear  ,  Cg  =  1.20  km/sec 

Bulk  ,  CQ  -  2.63  km/sec 


1 
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NICKEL 


Five  tests  were  performed  to  obtain  hugoniot  data  for  nickel. 
Three  of  the  tests  employed  nickel  impactors  and  two  impactors 
were  of  Fansteel-77.  The  nickel  was  purchased  under  specifi¬ 
cation  ASTM-B-160-61  as  99.5%  purity.  The  chemical  and  mech¬ 
anical  properties  of  the  material  are  listed  in  Table  11. 

The  experimental  data  are  tabulated  in  Table  12  and  displayed 
in  Figures  28,  29,  30  and  31.  The  least  squares  linear  fit  to 
the  data  is: 

U  =  4.456  +  1.555  u^  km/sec 
s  p 

Root  Mean  Square  Deviation  =  ±  0.012  km/sec  for  5  data  points. 

The  data  of  McQueen  and  Marsh (2) ,  Walsh (5) ,  and  Al’tshuler (18) 
are  also  displayed  in  the  figures  and  show  reasonable  agree¬ 
ment  with  the  present  data.  The  data  of  Al'tshuler  extends 
.over  a  wider  pressure  range  (1.1  to  9.2  megabars)  than  the 
General  Motors  data.  The  Al'tshuler  data  is  best  fit  by  a 
quadratic  curve: 

U  =  4.370  +  1.775  u  -  0.047  u  2  km/sec 
s  P  P 

Sigma  U  =  0.008  km/sec  for  4  data  points, 
s 
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TABLE  11  ■ 


CHEMICAL  AND  PHYSICAL  PROPERTIES  OF  NICKEL 


Chemical  Properties 


E lement 


Carbon 

Manganese 

Iron 

Sulphur 

Silicon 

Copper 

Nickel 


0.11 

0.26 

0.09 

0.005 

0.02 

0.02 

99.47 


Physical  Properties 

Yield  Strength 

at  0.2%  Elongation 

Tensile  Strength 

Density 


82,500  psi 


89,000  psi 
8.864  gm/cnT 


Poisson's  Ratio: 
Acoustic  Velocities: 


Longitudinal,  C^ 
Shear  ,  C 

s 

Bulk  ,  C_ 


0.300 

5.76  km/sec 
3.08  km/sec 
4.53  km/sec 
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Figure  29  Shock  Velocity  vs  Particle  Velocity  for  Nickel 
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M/NUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-68-13 

TYPE  304  STAINLESS  STEEL 

The  composition  and  mechanic*  I  properties  of  Type  304  Stain- 
’c  s  Steel  are  presented  in  Table  13.  The  measured  density 

is  7.905  gm/cm^.  The  hugoniot  data  over  a  pressure  range  of 

\ 

0.8  tr  4.3  megabars  are  listed  in  Table  14  and  are  presented 
graphically  in  Figures  32,  33  and  34.  The  linear  hugoniot 
fit  is  given  by: 


Ug  =  4.722  +  1.441  up  km/sec 


Root  Mean  Square  Deviation  of  Ug  =  t  0.023  km/sec  for  4  data 
points . 


Also  displayed  are  the  data  from  the  U.  S.  Army  Ballistics 

(19) 

Research  Laboratories  for  comtarison  with  the  present 
results.  Although  the  pressure  ranges  tested  are  barely  over¬ 
lapping,  the  extrapolation  of  the  present  dat  to  lower  pres~ 
sures  is  in  reasonable  agreement  with  the  Ballistics  Research 
Laboratories  results. 
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TABLE  13 

CHEMICAL  AND  PHYSICAL  PROPERTIES  OF  TYPE  304 
STAINLESS  STEEL 


Chemical  Composition 

Element 

Carbon 

Manganese 

Phosphorous 

Sulphur 

Silicon 

Nickel 

Chromium 

Molybdenum 

Copper 

Iron 

Qthwr  (CO) 


Weight  % 

0.065 
1.62 
0.029 
0.028 
0.49 
8.80 
18.  ?3 
0.14 
0. 17 
69.86 
0.070 


Mechanical  Properties 
Yield  Strength 
Tensile  Strength 
Hardness 
Density 


55,000  psi 
90,500  psi 
BHN  192 
7.905  gm/cro^ 


Poisson's  Ratio:  0.290 

Acoustic  Velocities:  Longitudinal,  C2  =  5.74  km/sec 

Shear  ,  Cg  =  3.12  km/sec 

Bu^  t  CQ  =  4.47  km/sec 


.  63 
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Figure  32  Shock  Velocity  vs  Particle  Velocity  for 
304  Stainless  Steel 
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Pressure  vs  Particle  Velocity  for 
304  Stainless  Steel 


Figure  33 
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TITANIUM 

Three  hugoniot  tests  were  performed  on  titanium  over  a  pressure 
range  of  0.4  to  2.7  Megabars.  The  chemical  and  physical  proper¬ 
ties  of  the  titanium  samples  are  listed  in  Table  15..  The  hug- 
oniot  data  are  presented  in  Table  16  and  in  Figures  35,  36,  37 
and  38.  The  linear  fit  to  the  hugoniot  is  given  by: 

U  =  4.692  +  1.126  u  km/sec 
s  p 

Root  Mean  Square  Deviation  of  U  =  ±  0.014  km/sec  for  3  data 

s 

points . 

In  the  figures  are  graphical  comparisons  of  the  data  of 
Krupnikov(20) ,  Walsh^5),  and  LASl/12)  with  the  present 
work  - 

Krupnikov:  Ug  =  4.8$  +  1.11  up  km/sec  (0.8  to  2.8  Mb) 

Walsh:  U  =  4.590  +  1.259  u  km/sec  (0.7  to  1.4  Mb) 

s  p 

LASL:  U  =  4.877  +  1.049  u  km/sec  (0.8  to  1.1  Mb) 

s  p 

A  surprisingly  large  spread  in  slope  and  intercept  is  seen, 
consolidation  of  all  the  data  above  yields  a  linear  fit: 


U  =  4.695  +  1.146  u 
s  P 

Root  Mean  Square  Deviation  of  O  =  ±  0.045  km/sec  for  13  data 
points,  deviating  but  very  little  from  the  General  Motors 
data  alone. 
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TABLE  15  1 


CHEMICAL  AMD  PHYSICAL  PROPERTIES  OF  TITANIUM 


Chemical  Composition 

Element 

Carbon 

Nitrogen 

Iron 

Oxygen 

Titanium 


0.02 

0.010 

0.25 

0.115 

99.60 


-  0.123 


Mechanical  Properties 


Yield  Strength 
at  0.2%  Elongation 


50,000  -  55,500  psi 


Tensile  Strength 
Density 


75,600  -  76,300  psi 

3 

4.508  gm/cm 


Poisson's  Ratio: 
Acoustic  Velocities: 


0.304 


Longitudinal,  C1  *  6.118  km/sec 


Shear 

Bulk 


,  C  *  3.246  km/sec 
s 

,  CQ  -  4.83  km/sec 
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*  BERYLLIUM 

Hugoniot  experiments  were  performed  on  two  types  of  beryllium, 
designated  as  S-200  and  1-400,  having  nominal  compositions  of 
2%  and  4%  BeO  respectively. 

The  chemical  and  physical  properties  of  the  beryllium  speci¬ 
mens  as  furnished  by  the  supplier,  Brush  Beryllium  Co.,  are 
listed  in  Table  17.  The  specimen  preparation  presented  some¬ 
what  of  a  problem  in  that  dust  and  small  fragments  of  the 
material  are  considered  toxic.  All  polishing  and  lapping  of 
samples  were  performed  in  closed  and  vented  work  areas. 

The  beryllium  presented  the  most  difficult  case  for  accurate 
hugoniot  measurements  with  a  restricted  target  size  because 
the  elastic  release  wave  velocity  is  relatively  high  and  the 
target  thickness  is  thus  necessarily  less  than  that  for  other 
materials.  (See  Figure  9  for  the  maximum  angle  of  intrusion 
of  the  side  rarefactions) . 

The  hugoniot  data  presented  in  Table  18  include  seven  tests 
on  S-200,  the  material  of  major  interest,  and  one  comparison 
test  on  1-400.  The  impactors  for  the  test  series  were  OFHC 
copper  (3  tests)  and  Fansteel-77  (5  tests).  The  data  are 
presented  graphically  in  Figures  39,  40  and  41.  A  linear 
fit  was  made  to  the  data  for  the  S-200  beryllium: 

U  *  8.390  +  0.975  u  km/sec 
s  p 

Root  Mean  Square  Deviation  of  Ug  =  ±  0.017  km/sec  for  7  data 
points . 
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TABLE  17 

CHEMICAL  AND  PHYSICAL  PROPERTIES  OP  S-200  AND 
1-400  BERYLLIUMS 


Chemical  Composition 


Material 

1-400 

S-^200 

Wt..  % 

Wt 

Beryllium  Oxide 

3.96 

2.00 

Carbon 

0.18 

0.126 

Iron 

0.16 

0.163 

Aluminum 

0.04 

0 . 054 

Magnesium 

0.02 

0.035 

Silicon 

0.04 

0.080 

Manganese 

0.01 

Beryllium 

95.77 

98.18 

Other 

0.10 

0.04  max 

Mechanical  Properties 

Yield  Strength 

56,000  psi  37,500  psi 

Ultimate  Tensile 

Strength 

66,700 

psi  57,700  psi 

Density 

1,881  gm/cirT 

1.857  gm/cm 

Poisson's  Ratio: 

0.055 

Acoustic  Velocities: 

Longitudinal,  C^ 

s 

12.83  12.916 

Shear  ,  Cg 

s 

8.80  8.86 

Bulk  ,  CQ 

= 

7.83  7.87 
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MATERIAL  -  0ERYLLIU4 


•  S-200  BERYLLIUM  GM  PQ  -  1.851 

♦  1-400  BERYLLIUM  GM  pQ  -  1.811 

0  LASL  (RZE.  12)  PQ  »  1.850 


LuiS- 

CL 

D 

in 

f^js. 

CL^ 

CL 


0-33  0*36  0-43  0-49  0-53  0-58 

SPECIFIC  VOLUME  -  CC/GM 


Figure  41  Pressure  vs  Specific  Volume  for 
Beryllium 
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The  1-400  beryllium  is  slightly  more  dense  1.3%)  them,  the 
S-200  beryllium  due  to  a  greater  percentage  of  the  heavier 
beryllium  oxide  compound.  The  data  point  is  slightly  above 
the  fit  of  the  S-200. 

A  comparison  of  the  data  with  that  of  LASL^2^  on  a  beryllium 

of  similar  BeO  content  shows  a  significant  difference  in  the 

slopes  of  two  linear  fits,  the  slope  LASL's  data  indicating 

a  somewhat  "stiffer"  hugoniot  than  the  present  data 

(U  =  7.998  +  1.124  u  ).  The  reason  for  this  difference  is 
s  p 

difficult  to  determine.  The  possibility  of  edge  rarefactions 
affecting  the  measurements  on  the  GM  specimens  was  checked  by 
firing  several  tests  with  specimens  substantially  thinner  than 
the  value  dictated  by  the  analysis  in  Section  II.  No  signifi¬ 
cant  differences  in  the  measured  values  were  noted. 
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AZ31B  MAGNESIUM 

Four  hugoniot  tests  were  performed  on  AZ31B  magnesium  tooling 
plate  (pQ  =  1.773  gm/cm3)  material.  The  nominal  composition 
of  the  material  is  Mg  96.0%,  Al  3.0%,  and  Zn  1.0%.  The  data 
shows  more  scatter  than  was  found  in  tests  with  other  metals. 

The  hugoniot  data  are  presented  in  Table  19  and  in  Figures 
42,  43  and  44.  The  hugoniot  is  described  by  the  linear 
equation: 

U  »  4.551  +  1.209  u  km/sec 
s  P 

Root  Mean  Square  Deviation  of  Ua  *  ±  0.083  km/sec  for  4  data 
points. 

The  figures  compare  graphically  the  measured  values  v?x  fch  the 

(21) 

AZ31B  magnesium  data  of  the  LRL  Compiler'  and  the  data  on 
99.5%  pure  magnesium  by  LASL^12* .  Linear  fits  to  the  ocher 
data  are : 

Ue  =  4.516  +  1.256  u  km/sec  (LASL,  p_  *  1.745  gm/cm3) 

5  p  v 

3 

U  =  4.648  +  1.198  u  km/sec  {LRL  Compiler,  p  =  1.78  gm/cm  ) 
s  p  o 
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Figure  44  Pressure  vs  Specific  Volume  for  Magnesium 
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PLEXIGLAS 

The  Plexiglas  tested  was  identified  as  Rohm  and  Haas  Type  II 
UVA.  Only  two  tests  were  performed,  which  served  to  provide 
a  check  on  the  extrapolation  of  the  data  of  Hauver'  taken 
above  the  phase  transition  at  ^  0.2  megabars,  (line  A  of  Fig¬ 
ure  45)  into  the  higher  pressure  region  investigated  by 
(23) 

Bakanova'  '  (line  B) .  The  data  are  listed  in  Table  20. 

Fits  to  the  above  data  are 

U  =  3.51  +  1.25  u  km/sec  (Hauver)  (2.8  <  u  <  3.5) 
s  p  p 

U  =  3.10  +  1.32  u  km/sec  (Bakanova)  (3.0  <  u  <  8.0) 

s  p  -  -  P 

Figures  45,  46  and  47  display  the  present  data  as  compared 
to  that  of  Hauver  and  Bakanova  with  the  line  representing 
Bakanova's  data  between  3.0  <  u  <  8.0  km/sec.  The  two 

p  t 

data  points  determined  in  the  present  work  are  in  good 
agreement  with  the  data  of  Bakanova  (difference  in  shock 
velocity  less  than  0.4%)  and  also  are  in  agreement  with 
the  extrapolation  of  Hauver' s  data  (difference  in  shock 
velocity  of  less  than  2%  at  a  pressure  of  1.0  megabar). 
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Figure  47  Pressure  vs  Specific  Volume  for 
Plexiglas 
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QUARTZ  PHENOLIC 

A  single  test  was  performed  to  locate  a  very  high  pressure 
point  on  the  hugoniot  of  quart 2  phenolic.  Lower  pressure  work 
was  performed  at  this  laboratory  under  Contract  AF04- (694) -807 
under  directorship  of  the  U.S.  Air  Force  Ballistics  Systems 
Division  and  is  reported  in  Reference  25.  All  of  the  data  are 
listed  in  Table  21  and  are  displayed  in  Figures  48,  49  and  50. 

A  phase  change  in  clearly  evident  beginning  at  a  pressure  of 
*  0.2  megabars.  The  data  point  obtained  is  on  the  higher  pres¬ 
sure  phase  of  the  hugoniot. 

The  linear  fit  to  the  upper  portion  of  the  hugoniot  is  given 

by : 

U  =  1.949  +  1  364  u  km/sec 
s  p 

Root  Mean  Square  Deviation  of  0  =  ±  0.016  km/sec  for  5  data 

points . 

The  initial  density  of  the  quartz  phenolic  is  1.80  gm/cm3. 

The  properties  of  the  high  pressure  phase  suggests  that  this 
may  be  related  to  the  properties  of  quartz  under  high  pressure, 
i.e.,  the  high  pressure  polymorph-stishovite '  '  .  The  tests 

were  made  with  layup  directions  of  the  quartz  cloth  in  the 
phenolic  matrix  both  parallel  and  perpendicular  to  the  shock 
front,  but  orientation  effects  on  the  wave  velocity  are  ap¬ 
parently  negiigable  at  these  pressures. 
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SECTION  V 
SUMMARY 


The  research  program  has  resulted  in  the  relatively  high  pre¬ 
cision  determination  of  the  high  pressure  hugoniot  equations 
of  state  of  eleven  materials.  The  estimated  error  for  shock 
velocity  measurements  with  the  experimental  system  was  0.2% 
to  0.6%.  Estimated  error  in  calculation  particle  velocity 
was  0.05%  for  like-like  impact,  but  somewhat  higher  for 
dissimilar  materials  impact  (see  Appendix  B) . 

Table  22  summarizes  the  results  of  the  experimental  program  in 

terms  of  the  constants  of  the  linear  hugoniot  equation 

U  =  C  +  Su  .  Also  included  are  the  pressure  ranges  mea- 
s  o  p 

sured  and  the  RMS  deviations  of  the  measurements.  Figures 
51  and  52  show  the  data  points  obtained,  displayed  in  the 
P  vs  ^  plane. 

Comparisons  have  been  made  of  the  data  with  that  of  other 
researchers.  Space  does  not  permit  the  inclusion  of  all 
data  available  for  comparison  and  the  data  quoted  is  selected 
principally  on  the  basis  of  pressure  ranges  studied  and  simi¬ 
larity  of  the  materials  in  terms  of  density  and  composition. 

The  onset  of  melting  in  the  shock  front  is  evidenced  by  only 
very  small  changes  in  the  measured  parameters,  shock  velocity 
and  particle  velocity.  Accordingly,  very  accurate  data  are 
required  to  substantiate  theoretical  predictions  of  the  melting 
phenomenon.  Although  the  linear  fits  to  the  high  pressure  data 
reported  here  are  useful  for  calculating  shock  wave  propagation, 
it  is  likely  that  several  of  the  materials  experience  melting 
at  the  higher  pressures ,  and  the  linear  approximations  become 
less  accurate. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  research  program  has,  by  the  attention  to  detail  necessary 
to  produce  high  precision  measurements  in  multimegabar  hugoniot 
determination,  provided  the  opportunity  to  work  out  the  diffi¬ 
culties  and  reduce  the  errors  in  a  relatively  new  experimental 
system,  the  use  of  a  light-gas  gun  for  shock  wave  measurements  - 
Although  some  problems  still  remain,  the  light-gas  gun  system  is 
now  better  developed  and  more  efficient  than  when  the  program 
began . 

On  the  whole,  good  agreement  was  found  with  the  data  of  other 
investigators.  However,  several  interpretive  problems  are 
still  unresolved.  The  behavior  of  the  aluminum  alloy  2024-T4 
in  the  pressure  range  1  megabar  should  be  investigated  fur¬ 
ther.  The  small  differences  in  the  slopes  of  aluminum  and 
copper  as  measured  by  gun-launched  impactor  and  by  explosive 
techniques  may  be  significant  and  need  further  investigation. 

The  large  differences  in  slopes:  .'f  the  beryllium  and  the  nickel 
hugoniots  when  compared  with  other  researchers  are  beyond  the 
estimated  experimental  errors  and  need  further  systematic  tests. 
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Numerous  helpful  suggestions  on  material  selection  and  data 
interpretation  were  provided  by  Dr.  S.  M.  Taylor,  the  contract 
technical  monitor,  and  by  Mr.  G.  E.  Hauver,  both  of  the  U.S. 

Army  Ballistic  Research  Laboratories. 
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INTRODUCTION 


The  fo'.lowing  is  a  detailed  description  of  the  data  analysis 
method  employed  for  the  determination  of  hugoniot  state  data 
derived-  from  the  experimental  system  described  in  the  main 
body  of  this  report.  The  first  part  of  this  appendix  pre¬ 
sents  the  equations  and  assumptions  employed  in  che  analysis 
and  the  second  part  describes  a  computer  routine  written  to 
perform  the  calculations. 


Experiment  Analysis 


Figure  A-l  illustrates  the  system  schematically  with  the  follow¬ 
ing  nomenclature: 


impact  velocity 

cm/p sec 

pin  cap  thickness 

cm 

pin  gap 

cm 

target  thickness 

cm 

pin  position  on  line  1-A 

cm 

impactor  tilt 

radians 

shock  velocity  in  target 

cm/ysec 
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t  =  <VhBtanv  ®>  tan  3  +  hB  bB  gB 

B  - - -  —  sec(^Q)  +  - -  + 


2  u 

spt  pt 


{A- 3) 


t  =  (X3"hCtan|  0)  tan  0  +  hC 


D  seo<v®)  +  u 


2  U 
spt  pt 


(A- 4) 


where  subscripts  1,  A,  B,  C  refer  to  pin  locations,  and 

Uspin  “  caP  s^oc^  velocity  induced  by  impactor 

Upin  =  pin  cap  partical  velocity  induced  by  impactor 

uspt  =  P^n  caP  s^oc^  velocity  induced  by  target 

upt  =  pin  cap  particle  velocity  induced  by  target 

Uspin  and  upin  are  evaluated  hy  the  impedance  match  solution 
for  the  direct  impact  of  the  projectile  material  into  the  pin 
material.  The  particle  velocity  derived  from  this  analysis  is: 


u 


_  B 


4AC 


2A 


(A-5) 


where  A  =  p  S  -  p  ,  S  . 

op  p  Kopm  pin 


(A-6) 


B=pC  +p.C.+2pCv 
*op  op  *opin  pin  op  op 


c  -  <>opv(Cop  +  Spv) 


and  U  .  =  C  .  +  S  .  u  . 

spin  pin  pin  pin 


>A~7) 
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the  impactor  hugoniot  being  given  by: 


C  + 
op 


S  u 
P  P 


(A-8) 


For  the  pin  gap  velocity  correction  of  the  wave  interaction  be¬ 
tween  the  specimen  and  the  rear  surface  pins  on  the  impedance 
match  solution  is  employed  again.  Although  the  equations  are 
similar  to  those  already  derived,  the  hugoniot  state  of  the 
specimen  must  first  be  estimated  and  then  determined  by 
iteration* 


The  first  approximation  of  the  shock  velocity  is  found  by  the 
following  procedure:  Subtracting  Equation  (A-6)  from  Equation 
(A-7)  we  obtain 


tan  0  = 


3  -  v 

x3  [ 


Vi1 


Vbl  .  ga'Tj 

Uspin  ^  upin 


A-9 


Neglecting  several  small  corrections,  the  approximate  shock 
velocity  in  the  target  at  Station  B  is  given  by 


°B  hB 


/  ^2  ^1 

3/  (tB~tl)  “  V~  tan  0  +  U  ~  +  Iu~ 
/  spin  pin 


(A-10) 


The  particle  velocity,  u' ,  at  this  shock  velocity  is  given  by 
Equation  (A-5) ,  vhere; 


A  =  p  S 
op  p 


B  =  p  C  +  p_, D '  +  2p  S  V 
op  op  ot  op  p 


(A-ll) 


C  -  p  vie  +  S  v) 
op  op  P 
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The  first  approximation  of  the  shock  pressure  is  given  by: 

P*  =  potD'u'  (A- 12) 

The  approximate  shock  and  particle  velocities  thus  obtained 
are  used  in  the  more  exact  solution  for  shock  velocity  given 
by  Equation  (A-13) .  The  results  may  be  iterated  several  times 
until  the  desired  convergence  is  reached.  With  the  convergence 
test  set  at  0.01%  usually  only  two  iterations  are  required. 


U  =  h„sec 
SB  B 


/ 

V 

/  f  tan  0 

X2  -hBtan  [l 

■ 

f  (t  -t,) - 

L_Z= 

V 


bl  gB 

—  + 


B 


U 


spin 


2  u  .  U 
pin  spt 


!b _ ; 

2  u  . 

ptj 


(A- 13) 


The  particle  velocity,  u  is  calculated  by  Equation  (k-S)  with; 

XT 


A  =  p  S 
OP  P 


B=pC  +p.U  +2p  Sv 
op  op  ot  s  op  p 


(A-14) 


C  "  popv(Cop  +  Spv) 


Pressure  is  calculated  from 


P  =  p  .  u  u 
ot  s  p 


(A-15) 


where  it  is  assumed  that  P„  =  0 

o 
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The  procedure  is  then  applied  to  analyze  the  data  from  pins  1, 

A  and  C,  A,  1,  and  B,  and  A,  1  and  C  using  the  appropriate  values 
of  g,  h,  b,  and  X. 

The  procedure  described  above  is  a  refinement  upon  the  simplest 
form  of  average  velocity  determination  (i.e.  distance/time) 
where  the  time  term  is  corrected  for  impactor  tilt.  The  above 
refinement  results  in  an  average  effect  on  Ug  of  0.2%  compared 
with  the  simplest  form.  The  overall  effect  on  a  shot  basis  has 
been  found  for  metals  such  as  copper  and  Fansteel-77  to  reduce 
data  scatter  about  a  linear  form  of  the  hugoniot  in  U3-up. 


Data  Analysis  Program  -  SHOVEL 

The  computer  program  described  below  is  written  in  Mark  I 
Fortran  for  use  on  the  General  Electric  Timesharing  Computer 
Service.  This  Fortran  is  very  similar  to  Fortran  IV  used  on 
many  computers  and  provides  some  convenience  for  adapting  to 
teletype  I/O. 

The  program  has  two  sections  and  one  subroutine.  The  first 
portion  of  the  program  is  merely  an  operator  convenience  to 
convert  all  input  data  into  the  required  units  and  formats. 
The  latter  portion  calculates  the  shock  velocity  at  each  sta¬ 
tion.  The  subroutine  COMPUT  is  in  the  form  of  the  quadratic 
solution  that  results  in  a  meaningful  value  of  u^. 
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Definition  of  Terms 


The  terms  presented  below  are  in  the  following  order: 

1.  Values  required  as  input 

2.  Intermediate  terms 

3 .  Output  values 

1.  INPUT  TERMS 


TERM  NAME  UNITS 


Shot  No. 

SHOT 

I 

D 

RUN 

Tr..pactor  Material 

MATLP (2) 

H 

Target  Material 

MATLT ( 2 ) 

H 

Impactor  C_ 

CZERO 

km/sec 

Impactor  S° 

SP 

Impactor  pp 

RP 

gm/cm^ 

Target  p  ° 

Pin  1-B  time 

RT 

gm/crn 

DTB 

ysec 

Pin  1-C 

DTC 

ysec 

Pin  1  height 

HI 

in. 

Pin  A  height 

HA 

in. 

Pin  B  height 

HB 

in. 

Pin  C  height 

HC 

in. 

Pin  1  gap 

G1 

in. 

Pin  A  gap 

GA 

in. 

Pin  B  gap 

GB 

in. 

Pin  C  gap 

GC 

in. 

Pin  1-C  Spacing 

XI 

in. 

Pin  1-B  Spacing 

X2 

in. 

Pin  1-A  Spacing 

X3 

in. 

D 

H2 

D 

H3 

Polaroid  target  fiducial  spacing 

M2 

cts . 

Impactor  to  Window  #1  Ref.  Spacing 

LI 

cts . 

Impactor  to  Window  #2  Ref.  Spacing 

L2 

cts . 

X-ray  Flash  Time  Interval 

TSI1 

nsec. 

X-ray  #1  to  Pin  C  Time  Interval 

TSI2 

ysec. 

Front  Surface  Shim  Thickness  (If  Used) 

SHIM 

in. 

Tilt  Time  Interval 

TILT 

ysec. 

Approximate  Slope  of  Target  Hugoniot 
(S=1.5  if  zero) 

S 

Pin  Combination 

COMBO 

H 

Computer  Run  No. 

RUNS 

I 

I  -  Integer,  D  -  Dummy,  H  -  Holorith  Field 
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2 .  COMPUTED  TERMS 


TERM  NAME  UNITS 


X-ray  1  Scale  Factor 

FI 

in/ct 

X-ray  2  Scale  Factor 

F2 

in/ct 

Window  Fiducial  Spacing 

D 

in 

Impact  Velocity 

W 

ft/sec 

Impact  Velocity 

W1 

ft/sec 

Impact  Velocity 

V 

cm/ysec 

D 

USTAR 

2 

Pin  and  Po 

RPIN 

gm/cm 

Pin  Material  Hugoniot  Slope 

SPIN 

Pin  Material  Co 

C0PIN 

km/sec 

Pin  Cap  Thickness 

Bl 

in3 

Specific  Volume  (initial) 

SVOL0 

cm  /gm 

Quadratic  Solution  Term 

A 

Quadratic  Solution  Term 

B 

Quadratic  Solution  Term 

C 

Particle  Velocity  of  Pin  Cap 

UPIN 

cm/ysec 

Shock  Velocity  of  Pin  Cap 

USPIN 

cm/ysec 

Impactor  Tilt 

THETA 

radians 

Firsv.  Approximation  Us  (B  Station) 

USB1 

cm/ysec 

First  Approximation  Up 

UBl 

cm/ysec 

First  Approximation  P 

PBl 

Mb 

First  Approximation 

RT1 

gm/cnr 

First  Approximation  Sound  Velocity 

C2 

cm/ysec 

Particle  Velocity  Induced  by  Target 

UPT 

cm/ysec 

Shock  Velocity  Induced  by  Target 

USPT 

cm/ysec 

Shock  Wave  Tilt  in  Target 

THETI 

radians 

3 

I 


<3 


S'. 

*3; 

ws 

v*; 


| 

1 


3 .  OUTPUT  TERMS 


B  Shock  Velocity 

USB 

km/sec 

B  Particle  Velocj  ty 

UB 

km/sec 

B  Pressure 

PB 

Mb 

B  Relative  Volume 

VOL 

o 

Specific  Volume 

SVOL 

cm  /gm 
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calculations  with  the  output  being: 


C  Shock  Velocity 
C  Particle  Velocity 
C  Pressure 
C  Relative  Volume 
C  Specific  Volume 


the  "C" 

station 

use 

km/ sec 

UC 

km/sec 

PC 

Mb 

V0L1 

3 

SV0L1 

cm  /gm 

The  "D"  station  calculation  employs  the  C  station  approximate 
velocity  and  corrects  for  the  time  interval  difference  between 
B  and  C.  The  output  terms  are: 


Shock  Velocity 
Particle  Velocity 
Pressure 
Relative  Volume 
Specific  Volume 


US 

km/sec 

U 

km/sec 

P 

Mb 

V0L2 

SV0L2 

cm  /gm 

Additional  output  terms: 


D  Station  Calculation  of  Shock  Tilt  THET3  Radians 

Calculated  Irapactor  Tilt  Time  Interval  TILT1  tsec 

%  Difference  of  (B-C)/C  Shock  Velocities  DIFF  % 


The  Mark  II  Time  Sharing  System  allows  chaining  of  programs  and 
SHOVEL  is  continued  as  SHOVEN  and  an  external  file  INFORM  is  used 
as  the  data  source.  Up  to  seven  files  may  be  linked  and  called 
in  process. 


*  *  v  :  ^  DEVELOPMENT  •  GENERAL  MOTORS  lOSPCHATIOtl 


MSL-68-13 


SHOVEL 


100  COMMON  U.  A#  B»C 

1  10  REAL  DIMENSION  MATLP<2>.  MATLTC2) 

120  INTEGER  SHOT. RJN. COMBO. RUNS 
1  30C 

1  40  200  READ£  1  *  40 > SHOT* RUN. MAIL  PC  1  >.NATLP<2>»MATLT<  l>.MATLTC2> 

1 50  REAO< I ) CZERO. SP. RP# RT. DT8. DTC.HI. HA. HB. HC.  G 1 »  GA» GB. GC. X 1 . X2. X3. 

1 60  AH2.H3.M2.L1.L2.TSI 1. TSI2. SKIM. TILT.S 
170C 

180  40  FORMATCSX* IS. 12. 4A6) 

j  900 

200  PRINT. "ENTER  PIN  COMBINATION1  C!  START  *  I*  A  START  •  A)  AND  RUN  NO." 
210  INPUT  30.  COMBO. RUNS 

g  Q0Q 

230  30  F0RMATC2AI) 

2  40C 

250CMPTV  COMPUTE  TILT  AND  VELOCITY 

252  PRINT  6 1 . CZERO. SP. RT. RP. DT8* DTC.HB. HC* GI. GA.GB. GC. XI . X2* X3. 

254  aTILT.S 

256  61  FORMATt////6FS.3.2F7.5.4F6.4//5F6.4///> 

260  c 

270  IFCLS ) 46. 48. 46 

280  46  Fla. 4756E-034F2*.4771£-034D*1 1*9951 

290  IF(L2>41. 42.41 

300  41  D1*<LI*F1+L2*F2aO)/I2.0 

318  VEL1*£O1/TSI1>*1.0EA06 

320  42  02* CL 1 *Fl ♦D*M2*F2*»  9537aSHIM)/ 1 2» 0 

330  VEL2*CD2/<TSI2-DTC>>**.0£*06 

3  40  IFCL2143.44.43 
3  50  44  GO  TO  45 

360  43  W*VELl4  W!*VEL2 

3  70  GO  TO  49 

380  45  W*VEL2f Wi*O.0 

390  GO  TO  49 

400  48  REA1KDW 

410  W1*>0.0 

4200 

430CTC  CONVERT  INPUT  TO  CGS 
440C 

450  49  V3< M/3.281 >*t*  0E+02 

460  CZERO* CZERO* 1 • 0E+054 USTAR* USTAR* 1 . 0E  >05 
470  DTB»DT8*1 • 0E«064DTC*DTC*I • 0E*06l TILT* TILT* 1 »0E**06 
n80  H1*H1 *2. 54IHA*HA*2. 54JHB*HB*2. 541 IHC*HC*2»  54 
490  RPIN  «8.41lSPlN«1.42lC0PIN*3.8E+05lBI*5.08E-03 
5  00  G 1  =  G 1 *2. 54# GA*GA*2« 544GB* GB*2*  544 GC«GC*2. 54 
510  X 1 * XI *2. 544  X2*X2*2*  544  X3* X3*2» 54 
529  THET3*0.0ISVOL0*1.0/RT4SVOL»0.04SVOL1*0.0 
S30C  ~ 

540CSV  COMPUTE  SHOCK  VELOCITY 
550C 

569  IF  CDTBJU.  It.  1000 
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1**0  Jj  —  u  O  '  1  , 


SHOVEL  CONTINUED 


5  70  11  UB*  0*  0<  PB*0«  0J  US3- 0»  01  VOL*  0»  0 
5  30  GO  TO  5 

5  90C 

600  1000  B=RP*CZ£RG*RPIN*C0PIN+2*RP*SP*V 
610  A=RP*SP-RPIN*SPINl C*RP*V*CCZERO*SP*V> 

620  CALL  COMPUT 

630  UPIM*UI USPIN*C0PIN*SPIM*UPIM 

640  THETA*<V/X3)*CTILT-<GA-G1  >/(2*UP!N>  > 

650  USB1=CHB-H1>/CDTB-X2*THETA/V> 

660  1001  B*RP*CZER0+RT*US81+2*RP*SF*VJA*RP*SP 

6  70  C*RP*V*<CZ£RO-»-SP*V) 

630  CALL  COMPUT 

690  U81  =  UJ  PB1*RT*U$B1*UB1 
700  RT1*<RT*US81 )/(US81-U81 > 

7 10  C2* USB 1*S9RT<. 49* (< USB 1-U8t )/USSl >**2> 

720  I FC  S) 50»  50>  51 

7  30  50  S*  1. 5 

7  40  51  B*RP!N*C0PI.N*RT1*C2*2*RT1*S*UB1*  A*S*RT1-SPIN*RPIN 
750  C*S*RT! *U81 **2*RT1 *C2*UBt+PBl 
760  CALL  COMPUT 

7  70  UPT*UJ  USPT-C0PIN  +SPIN*UPT 
780  THET1*  < USB 1/V>* THETA 

790  DELT1  =  DTB-THETA*<  X2-  CH8-H1  >*THET1  »/V 

800  DELT1=DELT1+B1/USPIN*G1/<2*UPIN>“BI/USPT-GB/(2*UPT> 

810  USB*<HB-H1 > *SQRT( THET1 **2* 1 7/DELT1 

820  3  IFC • 0001- < ABS<  < USB- US81 )/USB) > ) 1 002»  1 003» 1 003 

830  1002  USB 1* USB 

8  40  GO  TO  1001 

8  50  1003  B=RP*CZERO-*-RT*USB+2*RP*SP*V»A*RP*SPlC*RP*V*<CZERO*SP*V> 

8  60  CALL  COMPUT 

8  70  U9*U1PB=RT*USB*UBIV0L*1-U8/USB 
880  SVOL0*1*0/RT>SVOL=VOL/RT 

3  90C 

9  00  5  IF<DTC>6, 6.7 

9  10  6  USC*0.0JUC=0.  01  PC*0»  01  VOL  1*0»  0*  TILT* TILT* 1 • 0E+06I  US-0»  0; U*0»  0 

920  Ti L  T 1  * 0*  0 »  USB*  USB* 1 . 0E- 05  J  UB*  UB* 1  *  0£- 05 »  PB*  P8* I . 0Z- 1 2 

9  30  P*0.0lVOL2=0.0 

940  SVOLl*0.  0J  SVOL2=0»01DIFF*0.  0 

9  50  GO  TO  2 

9  60C 

970  7  8* RP*CZERO^RPIN*C0PIN*2*RP*SP*V 
9  30  A*  RP*SP“ RPIN*SPIN| C»RP*V* (CZERO  +  SP* V> 

990  CALL  COMPUT 

1  000  UPIN*UJUSPIN*C0PIN+SPIN*UPIN 
1  010  THETA* ?V/X3>*< TILT- <GA~G l  >/C2*UPI.N>  > 

1  020  USC1  =  (HC-H1>/(DTC-X1*THETA/V> 

1030  1004  B=RP*CZER0+RT*USC1+2*RP*SP*V*A*RP*SP 
1040  C*RP*V*<CZERO+SP*VJ 
1 050  CALL  COMPUT 
1060  UC1=U* PCI=RT*USC1*UC1 
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Shovel  CONTINUED 


J  873 
1  080 
1  090 
1  100 
1110 
1  123 
1  130 
I  140 
1  150 
1  160 
1  170 
1  180 
1  190 
1200 
1210 
1220 
I  230 
1  240 
1250 
1  260 
1  270 
I  230 
l  290C 
1  300C 
1310 
1  320 
1  330 
1  340 
1  350 
1360 
1  370 
1  330 
1  390 
1  400 
1  410 
1  420 
1  430 
l  443 
1  453 
1  460 
t  470 
1  430 
l  490 
1  500 
1510 
1  520 
1  530 
1  540 
1  553 
1  560 


R  T2- <  RT*USC l )/ ( USCS- UO! ) 

C3=*USC  I  *S9RT<  ♦  49*  C  (USCl-UCt  )/USCl  )**2> 

B=»  RPIN*C0P1N+RT2*C3*2*RT2*S*L‘C1  J A*S*RT2-  SPIN*RPIN 
S*RT2*UC I **2*RT2*C3*UC 1*PC 1 
CALL  COMPUT 

UPT1*U*  USPT1*C0PIN*SPIN*UPT1 

THET2*<USC1/V>* THETA 

DEL  T2*  DTC- THETA*  <  X 1 - <  HC- H I  »  *  THET2 > / V 

DELT2aDELT2+Bl /USPIM+61/ t2*UPIM)»81/USPTl“QC/C2*UPTt > 

USC» ( HC- H l >  *  S9RT< THE  T2**2+ 1 ) /DEL T2 
4  IF<  • 0001- < ABSC  ( USC- USC 1 >/USC> ) ) 1 005* 1006*1006 

1005  USC 1 * USC 
GO  TO  1004 

1 006  B*RP*CZER0+RT*USC*2*RP*5P*Vl A«RP*SP*C*RP*V*<CZERO+SP* V> 
CALL  COM PUT 

UC*Ul PC*RT*USC*UCl  V0L1» 1-UC/USC 

UB«UB*  1  •  0E-5*  USB* USB*  1 . 0E-  05 1  PB«  PB*  I  .  0E-  1 2 

uc-  UC*  1 . 0E- 051  USC«  USC*  I  *  0E-  05  J  PC*  PC*  1  •  0E-  i  8 

DEL  T3 <  0SLT1 - DEL  T2 ) * 1 . 0£*  06 

T!LT3TILT*1.0E*06 

TILT1*<DTB-0TC»*X3/<X2-X1 >*1.0E*06 
DIFF«0.0*SVOL1*VOL1/RT 


IF<DTB>70* 70,21 

70  TILTl30.0JUa3«  01  USa0*  01  P»0*  0t  VOL2*0*  0 

GO  TO  2 

21  US1*<HC-H1)/0TC 
BaRP*CZER0*RT*USl+2*RP*SP*VlA*RP*SP 
C*RP*V*CC2ER0*SP*V) 

CALL  COMPUT 

UI*UlPl3RT*USl*Ul 

1007  RT3*<RT*USl >/<USl-Ul > 

C43US1 *S9RT( »49*CtUSt-Ul >/USl>**2) 

IF<S>22* 22*23 

22  S*  I  •  5 

23  83RPIN*C0P[M+RT3*C4*2*RT3*S*U1IA*S*RT3-SPIN*RPIN 
CaS*RT3*Ul **2*RT3*C4*Ul*PI 

CALL  COMPUT 

UPT*Ut  USPT“C0PIH-f SPIH*UPT 

0T1LT=DTB-<HB-H1 >/US1*81/USPIN+G1/(2*UP!H)-B1/U3PT-GB/<2*UPT) 
THET33V*DTILTAX2 

SHTIM3DTC-  THET3*<X1-(HC-Hl  >*USl*TH£T3/V)/V+Bl/WSP!N+Ql/<2*UPIiM> 
>-Bl/USPT-GC/(2.*UPT) 

US3 (HC-H 1  >  *S9RTC  1  •  0*  C  US  l * YHET3/ V) **2) /SHTIM 
83RP*CZ£R0*RT*US+2*RP*SP*Vl A*RP*SP 
C=  RP* V*( CZERO*  SP*V> 

CALL  COM PUT 
U3Ul  P=*RT*US*U 

l F< .0001- CABS<<US-US1 >/US>  >> 1008* 1009* 1009 


oivau^MiH  i  *  general  motors  corporation 

MSL-68-13 


SHOVEL  CONTINUED 


1570  1008  USl*USlUl*UlPl*P 
1 530  60  101007 

1590  1009  VOLasl-U/(JSSU*U*1.0e-0S!US«US*1.0E-05lP*P*'»0£-lS 
1600  SV0L1»V0LI/RT*  SV0L2*V0L2/RT 
1  610  0I?T«<  ABS< < USB-USJ/USG)  >*  1 00*  0 
1620  SUSE  5H0VEN 


S WOVEN 


1700C  OUTPUT  STATEMENTS  FOR  SHOVEL 
S  710C 

1720  2  PRINT  60#  S40  T»  RUNS#  COMBO# MATL TC  I  >»MATLTC2>#  MATLPC  I  >»  MATLPC  2># 

»  730  ♦HT#RP#SVOL0#W#TH£T3#V#TILT»W1#TILT1.DIFF# 

*  740  +USB#  UB»  PB»  VOL#  SVOL#  USC# UC#  PC#  VOL  t  #  SVOL 1 »  US#  LI#  P#  V0L2#  SV0L2 
1  750C 

2760C  OUTPUT  FORMAT 
l  770C 

1  780  60  FORM AT< //////// 10X#  9HS40T  NO.  # 1 5# 1 0X# 8HRUN  NO.  #AI// 

1790  + 1 0X#  1  7HPIN  COMBINATION*  #A1#1X#7H-  START#// 

1800  +5X# 14WTARGET  MATL.-  #2A6# 

1810  +4X. 13HPR0J.  MATL.-  #2A6//8X#9HDEN5ITT-  # F6»  3# IX#  5HGM/CC#  8X» 

3  820  +9HDENSI TY-  #  F6«  3#  1 X# 5HGM/CC//8X# 8KV0LUME-  »F6*  4#  1 X#  5HCC/GM# 

1830  *///// 5X# 1 8HPR0JECTILE  VEL.  *  #Et2»5# 

1840  ♦ 1 X#  6HFT/SEC. 4X» 1 3HSH0CK  TILT  *  #F9»6//23X#EI2.5# IX#  6HCM/SEC# 

1850  *4X»  13MTILT  MEAS.  »  *F9.6//23X»E12»5»  IX#  6HFT/SEC#  4X# 

I860  +13HTILT  CALC*  =  #F9. 6///5X#21HSH0CX  VEL-  SCATTER  =  , F7. 3# 1 X# 1 H J# 
1870  +//////2X#  4HS'; A.  4X#  1  0HSHOCX  VEL.  i'X#10HPART.  VEL.  7X>  8H PRESSURE 
1  880  4-4X#  4HV/V05X#6HVOLUME//10X#6HKM/$EC#8X#6HKM/SEC»  1 1X#8HMEG^  "'RS# 
1890  + 1 3X#  5HCC/GM//// 

1  900  +2X#  1  HB»  6X»  F8.  4#  6X# FS*  4#  8X#  FI 0.  6j  1 X#  P9.  5#  1 X#  F8«  5// 

1  910  +  2X#  1HC#6X#F8.  4#  6X»F8.  4#  8X»  FI0.  6#  1 X.  F9.  5#  1X#FS.5// 

1920  +2X#  1HD#6X#F8.  4#  6X#  F8.  4#  8X#  F10.  6»  1  X.  F9. 5#  1 X#  F8#  5////// > 

1930  GO  TO  200 

1940  END 

1950  SUBROUTINE  COMPUT 

I960  COMMON  U#A#8#C 

1 970  U»<B-S9RTCB**2-4*A*C>>/<2*4> 

2  980  RETURN 
1990  END 

2000  SFILE  INFORM 
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APPENDIX  B 

ANALYSIS  OF  SYSTEMATIC  MEASUREMENT  ERRORS 


An  analysis  of  systematic  errors  in  the  experimental  system 
* 

'<•/  i'i  performed.  Each  of  the  22  parameters  involved  in  the  cal¬ 
culation  of  a  hugoniot  equation  of  state  point  from  information 
obtained  using  the  light  gas  gun  system  was  considered  and  an 
estimated  error  was  calculated.  These  errors  were  combined 
mathematically  to  yield  an  average  systematic  error  for  the 
hugoniot  parameters  calculated?  pressure,  particle  velocity, 
shock  velocity,  and  volume.  Errors  in  measurement  of  impact 
velocity  were  also  calculated. 

Two  cases  were  considered:  (1)  dissimilar  impact,  in  which  the 
hugoniot  point  is  calculated  from  the  impact  velocity,  the  shock' 
wave  velocity,  and  the  measured  hugoniot  for  the  impactor  and 
(2)  symmetric  impact,  in  which  only  the  impact  velocity  and  the 
shock  wave  velocity  are  necessary  for  the  calculation. 


Since  the  error  analysis  is  dependent  on  the  shock  impedances 
of  the  specimen  and  the  impactor  as  well  as  on  the  pressure 
attained,  a  complete  error  analysis  could  be  performed  for  each 
data  point.  This  formidable  task  was  avoided  by  calculating 
errors  involved  in  an  "average"  data  point,  and  briefly  examin¬ 
ing  the  errors  at  the  limits  of  the  system.  As  the  variance 
in  errors  at  the  limits  do  not  substantially  change  the  values 
for  the  average  case,  only  the  average  case  is  reported. 


The  parameters  operated  upon  involved  a  measurement  made  on  copper 
with  a  Fansteei-77  impactor  or  a  copper  impactor  launched  at  6km/sec. 
The  method  is  a  statistical  procedure  called  the  analysis  of  variance. 


The  authors  are  indebted  to  Mr.  W.  L.  Rearick  of  Systems 
Analysis,  Manufacturing  Development,  for  his  assistance  in 
the  formulation  of  the  system  error  analysis  procedure. 
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By  definition  the  standard  error  of  the  variable,  x,  is  given 
by: 

°SE  '  V  <j2xi''n  <a) 

where  ox^  is  the  deviation  of  the  ith  term  from  some  f (x) . 

When  ox,  a2,  an  are  mutually  dependent,  the  variance  of 

x  is  then: 


'xl 


+  a 


x2 


+  a 


xn 


(b) 


.  2 

with  denoting  the  nth  element  component  of  x. 

In  the  present  case,  a  a  (or  standard  deviation)  is  determined 
for  the  system  components  of  all  the  mutually  dependent  variables, 
listed  in  Table  1-B* 

The  variances  of  mutually  independent  variables  may  be  determined 
using  the  rules  presented  below.  By  definition  the  mean,  y  ,  of 
a  variable,  x,  is  given  by: 

Px  =  (X1  +  x2  +  •  •  •  xn>/n  (c) 

where  x^,  x2,  .  .  • /  xn  are  element  components  of  a  statistical 
distribution. 


If  the  variance  of  x  and  y  are  small  compared  their  means,  i.e., 
2 

o  <<  \i  ,  and  independent,  the  following  rules  may  be  employed 

X  X 

to  determine  the  variances  of  the  sums,  products  and  quotients 
of  variables  x  and  y: 
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Addition:  a2  =  0  2  +  a  2 

x±y  x  y 

Multiplication:  c  2  =  o  2a  2+c  2y  2+a  2p  2 

xy  x  y  x  *y  y  Hx 


Division: 


°xW^2 


(f) 


where  in  (f)  , 


y  >>0 
y 


The  present  analysis  was  carried  out  using  the  tabulated  values 
of  y  and  0  of  Table  1-B.  The  variance  of  each  term  of  equations 
1  through  4  was  determined  and  then  combined  to  yield  the  re¬ 
sultant  variance  for  each  operation. 


Impact  Velocity:  v  =  y/z 


(1) 


Impactor  Tilt: 


tan  6  =  -  (At, - ) 

x,  t  2-u  .  ' 
3  pin 


(2) 


Shock  Velocity:  U  =  h  sec  (=--)/[  (T. -T  )  - 

a_  «  v  I  n 


n 

x  — h  tan 9  tan  ( **■  8 ) 
n  n _  y  ' 

v 


U.  2*u  . 
spin  pin 


(3) 


'n 


2*upt 


] 


Symmetrical  Condition  Particle  Velocity: 


Up  =  V2v 


(4) 
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To  determine  the  variance  resulting  from  the  impedance  mismatch 
solution  for  the  particle  velocity ,  it  is  necessary  to  know  the 
variances  of  the  terms  of  the  hugoniot  equation  of  the  standard. 
Since  the  hugoniot  equation  is  itself  a  statistical  inference, 
the  following  approach  was  used: 


y  =  s?  J2i 

up  1  n 


where  n  is  the  number  of  data  employed  in  the  least  squares 
analysis.  The  variance  is  then  given  by: 

0  2  =  e”  (u  . -y  )2 

up  1  v  pi  Mup 

For  an  equation  in  the  form: 

U  =  C  +  S 
s  o  up 

the  variances  of  the  zeroth  and  first  order  terms  using 
equation  (a)  are  given  by: 


2  _  °se 

s  “  2 

P  0 
r  up 


The  analysis  of  variance  is  then  applied  to  equation  5a  through 


Impedance  mismatch  particle  velocity: 


u  = 
P 


B-Yb  -4AC 
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A  =  Pop  '  Sp 


E  =  p  C  +p  ,U  +2*p  -S  *v 
op  op  Kot  s  Fop  p 


C  =  p  C  v+p  S  v‘ 
op  op  op  p 


To  determine  the  variance  of  the  square  root  operation,  the 
binomial  expansion  is  used;  i.e. 


u2  +  a2  =  2^ 


The  variances  are  then  determined  for  the  remaining  hugoniot 
relations . 


Shock  Pressure: 


P  =  p_.  U  u 
ot  s  p 


Shock  Specific  Volume 


fe.  v  = 

ie .  V  y. 

s 


Table  2-B  presents  a  summary  of  the  systematic  error  analysis. 
The  percent  error  column  shows  the  marked  difference  in  the 
accuracy  with  which  measurements  can  be  performed  on  symmetric 
(Case  2)  versus  dissimilar  material  (Case  1)  experiments.  The 
larger  errors  involve^  in  the  dissimilar  material  tests  (Case  1) 
concern  the  errors  in  the  hugoniot  of  the  impactor  material. 

For  tests  using  impactor  materials  with  precisely  measured 
hugoniots  these  errors  will  be  correspondingly  smaller. 
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19.  ABSTRACT 


An  experimental  technique  is  utilized  in  which  a  light-gas  gun  is 
used  to  launch  flat  impactor  plates  to  high  velocities  8  km/sec) 
at  specimens  suspended  at  the  muzzle  of  the  gun.  Impact-induced 
shock  waves  at  pressures  to  5  megabars  are  recorded  and  are  used 
to  determine  the  shock  state  in  the  specimen.  The  ability  to  launch 
unshockeJ,  stress-free  flat  plates  over  a  wide  and  continuous  velocity 
range,  coupled  with  the  ability  to  launch  impactor  plates  of  the  same 
material  as  the  target,  results  in  hugoniot  measurements  of  relatively 
high  precision.  Measurements  were  made  on  Fansteel-77  (a  tungsten 
alloy),  aluminum  (2024-T4) ,  copper  (OFHC,  99.99%),  nickel  (99.95%), 
stainless  steel  (type  304),  titanium  (99.99%),  magnesium  (AZ31B) , 
beryllium  (S-200  and  1-400),  uranium  (depleted),  plexiglas,  and 
quartz  phenolic.  The  results  are  compared  with  those  of  other 
researchers . 

Deviation  from  linear  shock  velocity  -  particle  velocity  was  found 
in  aluminum  beginning  at  ^  1.0  megabars,  probably  attributable  to 
melting  in  the  shock  front. 
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